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DEVELOPMENT AND EVALUATION OF SCATTEROMETER 
DATA PROCESSING ALGORITHMS 

1.0 INTRODUCTION 

In early 1976 the NASA Johnson Space Center (JSC) Earth 
Observation Division contracted with the Remote Sensing Center 
to develop and implement the hardware and software required to 
process raw radar measurements acquired using the 0.4 GHz, 

1.6 GHz, and 13.3 GHz scatterometers to normalized radar cross 
section. These scatterometers had been taken out of retirement 
during 1976 and mounted on the NASA C-130 aircraft to be flown 
in support of the NASA funded Agriculture Soil Moisture Experi- 
ment (previously termed the Joint Soil Moisture Experiment) . 

A data processing system was developed and implemented by 
Texas A§M University in 1976. The processing system was 
structured to consist of two phases. First, digital magnetic 
tapes of raw scatterometer measurements and ADAS/NERDAS data 
were generated using the Remote Sensing Center (RSC) TI-980 
mini -computer, a four channel analog-to-digital converter con- 
structed by the RSC, an ADAS/NERDAS discriminator constructed 
by the RSC, and a fourteen track analog recorder supplied by 
NASA Johnson Space Center as government furnished equipment. 
These tapes were then processed to normalized radar cross 
section using a software package implemented on the TAMU 
Amdahl 470. This processing system was documented in the 
initial "final report" on the contract (1J. The program was 


unique in that it utilized Hilbert transforms to sign sense 
the scatterometer data and discrete digital Butterworth filters 
to select the component of the return correspondence to spe- 
cific incident angles. Detailed flow charts for this system 
are deluded in Appendix A. 

Radar scatterometer data acquired during the 1976 Joint 
Soil Moisture Experiment at the LACIE supersite in Finney 
County , Kansas and a test site at Lawrence, Kansas was pro- 
cessed using this processing system in 1977. In evaluating 
these data several questions were raised concerning its 
validity. As a result, considerable effort wa. c devoted to the 
evaluation of data acquisition, data pre-processing and final 
data processing. During this evaluation it was discovered 
that the analog tape recorder that had been provided as 
government furnished equipment by NASA Johnson Space Center 
(JSC) could provide only marginal performance and was inade- 
quate for processing scatterometer data. These problems were 
summarized in the quarterly progress review presented at 
NASA/ JSC on September 22 and 23, 1977. The pertinent portions 
of this presentation are contained in Appendix B. As a result, 
another more suitable analog tape recorder was provided to 
TAMU by NASA/JSC for processing scatterometer data. 

Also during this evaluation, several other sources of 
error were discovered. The most significant of these were: 

1) poor performance of the Hilbert transform at low incident 
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angles, 2) use of average flight parameters for each flight 
line to compute cell geometry and area, 3) poor and/or inade- 
quate flight log documentation leading to misinterpretation 
of sensor parameters, 4) poor communication between sensor 
system engineers at NASA and TAMU leading to the use of incor- 
rect and/or outdated sensor calibration constants, and 
5) inconsistent identification of the sign and cosine channels 
of each scatterometer, and misleading identification of 
channel assignments of these signals to tape recorder tracks 
on the analog magnetic tape labels. 

After the above mentioned problem areas were discovered 
and many of them resolved, NASA/ JSC began to implement a 
software processing system in hopes of being able to process 
future radar data in-house at less cost. TAMU was given 
additional funding by the NASA/ JSC Earth Observation Division 
to support the establishment of this software system by the 
Engineering and Development Division and evaluate its perfor- 
mance in conjunction with the continuing evaluation of the 
TAMU software system. TAMU was also given permission to expand 
its processing capability to include the newly acquired 
4.75 GHz scatterometer. Within this same time frame TAMU was 
contracted by NASA/JSC Engineering and Development Division 
to develop and design a real time hardware on-board processor 
for the 1.6 GHz and 4.75 GHz scatterometer systems. Therefore, 
the hardware processor design and software processor evaluation 
proceeded simultaneously. 
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The simultaneous activities on the hardware and software 
processing systems proved to be extremely beneficial. Pri- 
marily as a result of the new techniques that were being 
developed for the hardware processor [2], two new software 
systems were developed for processing scatterometer data. The 
Hilbert transform and discrete digital filter approach was 
abandoned in favor of approaches that relied on discrete 
Fourier transforms for sign sensing and filtering. The major 
advantage of this approach is that the filters can be modified 
during data processing as a function of aircraft parameters 
in order to maintain constant resolution and incident angle. 

One system developed by TAMU relies on real valued 
discrete Fourier transforms of each scatterometer quadrature 
channel. Sign sensing is accomplished by algebraic manipula- 
tions of each of these transforms. This system is identical 
in processing procedure to the system developed by NASA/ JSC. 

The other processing system developed by TAMU has the identical 
processing procedure that will be used in the on-board 
scatterometer hardware processors. This system relies on a 
single complex valued discrete Fourier transform to sign sense 
and filter the quadrature signals [2]. This system is the 
best in terms of speed of operation since 1.6 GHz, 4.75 GHz, 
and 13.3 GHz data only half of the number of Fourier transforms 
must be computed. And it provides identical results. 
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The evaluation phase of the effort involved much personal 
communication between TAMU and NASA/ JSC employees. The many 
performance tests and equation derivation verifications that 
were painstakingly executed will not be included in this report. 
Suffice it to say that agreement was finally achieved between 
the software processing system at NASA/JSC and the two pro- 
cessing systems developed at TAMU utilizing Fourier transform 
techniques. However, the programs developed at TAMU and 
NASA/ JSC are not identical in their manner of construction 
or execution. Parameters such as resolution, integration time, 
percent of data utilization, etc. are different for the two 
processing systems. Users of scatterometer data should be 
aware of these differences since they could have important 
implications in the analysis of certain data sets. 

In the following text the software processing systems 
developed by TAMU as a final product will be described. 
Differences in these programs and the NASA/JSC program will 
be pointed out where appropriate. It should be noted that the 
TAMU software system utilizing the single complex discrete 
Fourier transform for sign sensing and filtering operates in 
exactly the same manner as the TAMU design for the real time 
hardware processor [2]. 
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2.0 GENERAL SYSTEM DESCRIPTION 

The software system described in this document is the 
final product of the project and is currently being used 
exclusively for all radar scatterometer data processed at TAMU. 
The general processing procedures are described in Figures 1*4. 
The first step in processing scatterometer data is to digitize 
the raw radar measurements as well as the ADAS or NERDAS data 
strewn, Figure 1. Both quadrature channels of radar measure- 
ments are digitized simultaneously for one radar frequency 
at a time. ADAS or NERDAS data are discriminated and digitized 
at the same time and a computer compatible digital magnetic 
tape is generated that contains the radar measurements and 
time correlated ADAS/NERDAS data. The digitization rate of the 
A/D converter is software controlled and depends on the radar 
frequency being digitized. 

The percentage data utilization is a function of the A/D 
process and varies for each radar frequency. One hundred per 
cent utilization would mean that all of the raw radar measure- 
ments are digitized. This does not occur for either the TAMU 
processor or the NASA/JSC processor. However, this does not 
mean that the final output product does not provide continous 
ground coverage. 

The minimum integration time of the processor is set by 
the record lengths and A/D rate used during the A/D conver- 
sion process. Both the TAMU and NASA/ JSC processor use the 
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STEP 1: Data Digitization 


Raw Radar Data 




Raw Radar Data 
Digital 




1 Record/Cell 

4112 Bytes/Record (Includes 
aircraft parameters and 
radar data) 

Radar data - 1024 words/record 
for each channel (sine and 
cosine) 


ADAS/NERDAS: Time, 

altitude, velocity. 

pi tch , roll, drl ft 

RADAR Frequency 

A/D Rate 

Percentage Data Utilization 

13.3 GHz 

25 KHz 

16.4 

4.75 GHz 

10 KHz 

37.2 

1.6 GHz 

5 KHz 

51.2 

0.4 GHz 

2.5 KHz 

68.2 


Figure 1. Simplified flow diagram of the first step of 
the radar data processing procedure-Data 
Digitization. 









sane record lengths and A/D rates. The minimum integration 
time determines the smear distance of a resolution cell. 

Although the basic integration times are the same between 
NASA/ JSC and TAMU, NASA/ JSC defines resolution differently than 
TAMU thereby causing the actual ground coverage obtained from 
the two processing systems generally to be different. This 
is discussed later. 

Each radar frequency and polarization combination must be 
digitized individually. When this is complete , the 9 track 
computer compatible tape is used to input the digital data to 
the TAMU Amdahl 470 computer where the main processing soft- 
ware resides. The Amdahl 470 processes the measurements to 
normalized radar cross section using the program SCATTER. CZTA78 
and produces both punched card and line printer output (Figure 2). 
The processed data are then put into report format for distri- 
bution using the TI-P80 minicomputer, plotted for analysis 
purposes, or avevaged as desired (Figure 3). 

The principal component of the processing system is the 
program SCATTER.CZTA78. This program sign senses, filters, and 
computes the normalized radar cross section. A general flow 
chart of this program is shown in Figure 4. 

In evaluating the processing system, each step of the 
process was investigated individually. The accuracy of the 
digitization process was tested and the ADAS and NERDAS 


STEP 2: Sigma Zero Computation 




Data 

Cards 


Non-Al Igned 
Data List 


Figure 2. Simplified flow diagram of the second step 
of the processing procedure - Sigma Zero 
Computation. 
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Product Production 



Simplified flow diagram of the third step of the data processing 
procedure with indications of the type of output product available. 










Figure 4a). Simplified Flow Oiag 





of Program SCATTER. CZTA78 








Figure 4.b) Antenna gain, beamwldth and 

roll-off function modules are loaded 
Into main proqram Scatter at time of 
program execution for the frequency 
and polarization combination being 
processed. 
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discrimination tested. These were found to work satisfactorily 
and compare to NAS A/ JSC results for the same input data. The 
quality of the 9 track digital tape is tested each time a set 
of radar data are processed by executing a fast Fourier trans- 
form of the first few records and comparing the spectrum to 
that observed (and photographed) on the analog spectrum 
analyzer when the analog tapes are digitized. All of these 
functions proved to bo satisfactory after the original analog 
tape recorder problems were rectified. 

The major effort in the evaluation process proved to be 
the validation of the main program component, SCATTER. CZTA7 8. 
This program basically evaluates the radar equation as given by 
equation (1) for the normalized radar cross section. The 
assum tions upon which equation (1) is based were evaluated, 
and the techniques and accuracy of computing each of the terms 
in equation (1) were evaluated. The details of each of these 
evaluations will not be given, only a summary of several major 
points. 

3.0 GENERAL RADAR EQUATION SOLUTION 

At the outset of this report period two forms of the 
general radar equation were used to estimate the radar cross 

A 

sections, o , from the scatterometer measurements. At TAMU 
the form used is 

o-o = 62 ? L 2c S! d) 

A.** Kx <**(*£/% /t 
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where 


X * wavelength 
* cable loss 

K s ■ system constant, calibration 
Z(f) ■ roll-off function 
G r G t * antenna two-way gain 
P r * receiver power 
P t ■ transmitter power 
R = range to ground cell in view 
A = area of ground cell in view 

At NASA/JSC the form used was 


CT 6 - 


urfr c±_ V/SV _ \(£» i V A <- Vans, I 


( 2 ) 


where 


H 

V 


A 


c 


$ 


BW C 



* 

4> 


rms volts corresponding to the received band 
width 

rms volts corresponding to the calibration 
bandwidth 

aircraft altitude 

aircraft velocity 
A0 

COS (-J-) COSlJ>COS<j> 

antenna beamwidth 

bandwidth of the calibration signal 

instantaneous bandwidth corresponding to the 
i t " incident angle 

roll angle 

drift angle 
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The two expressions are related and equivalent if it is 
assumed that the ground cell area, A, in equation (1) is 
given by 

** " (3) 

where 

th 

8. * instantaneous viewing angle to the i ground 
1 cell 

A8 = the angular width corresponding to the BUL 

The above assumption which connects equation (2) to equation 
(1) is, unfortunately, valid only if the value of roll and 
drift remain equal to zero. This is almost never true, there- 
fore, the authors of the orginal NASA/ JSC system were forced 
to make another assumption to cover this situation; i.e., in 
the presence of aircraft roll and drift the area becomes 

(4) 

The validity of the above assumptions was extensively analyzed 
at TAMU [3]- [4]. It was found that equation (3) (first 
assumption) waa accurate to within ±1% for the case where 
H 8 0; however, for <J>^0, equation (4) (second assumption) 

yielded erroneous values of ground cell area which were off by 
as much as 7% for nominal aircraft perturbations. On the basis 
of these analyses it was concluded that the more general form 
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of the radar equation, expression (1) above, would be the best 
form to use for data processing. 

The results of this analysis also showed that the area 
calculation should be done using the expressions and techniques 
that have always been used at TAMU; i.e., 

h/'C y* -y>) ( 5 ) 

Ces & 

as defined in reference [4]. The NASA/JSC program was modified 
to compute area in the same manner for purposes of consistency 
between the two software processors and the hardware processors 
which are being implemented. 

4.0 DOPPLER FREQUENCY SHIFT CORRECTIONS 

4.1 Effect on Range and Area Computation 

4 

The value of R /A (equation (1) above) obviously makes a 

A 

substantial contribution to the value of o°; therefore, a 
careful evaluation of the technique for computing these values 
was done. It was found that adjustments had to be made to 
account for doppler filter center frequency and doppler filter 
band-edge shifts associated with the discrete nature of the 
Power Spectral Density (PSD) lines. For example, the center 
frequency of the filter corresponding to a particular viewing 
angle is given by 


where 0 L is the viewing angle with respect to the aircraft 

gravity vector. Now, as illustrated in Figure S, this value 

may or may not fall as the center of one of the PSD lines. 

Similarly, the bandwidth, BW^, is actually some integer set 

of PSD lines, and thus the filter band-edges do not actually 
BW, 

fall at fd+ — , but in fact fall at the edge of one of the 
PSD lines. The true bandwidth may be calculated as 


where the notation [ ] j means taking the integer value, and 
Af is the spectral line width given by Af » f /N where N is the 
number of lines in the PSD and f is the sample frequency used 

9 

by the A/D converter. The true center frequency is then 
given by 

+ ( 8 ) 

whenever the expression 



( 9 ) 


is odd, and by 

At, * *?[£&/*$ +*(L +*■€ 

whenever NF^ is even. 
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Figure 5. Effect of discrete nature of the DFT on filter 
center frequency and bandwidth. 




Note that at most, the center frequency or the band-edge 
will be off by A */2. The influence of this perturbation on 
the value of range, R, can be evaluated. Range is given by 

4V = At <“> 

Since and f d are related to equation (6) , we can write 





where R£ is the range to the repositioned ground cell with the 
center frequency off by 2 . It can be shown that this func- 
tion is maximum when the value of Xf„ is maximum. This maximum 
occurs at the lowest frequency, 0.4 GHz. Also, the largest 
percentage difference between Rj and R! occurs at the smallest 
0 L ; i.e., normally S degrees. Thus, the ratio in equation (6) 
can have values as large as 1.07. This implies that if the 
center frequency shifts are not accounted for in the calcula- 

4 

tion of R the error can be as much as 33% when processing 
0.4 GHz data. When this is applied to equation (1) to compute 

A 

a , the result could be off by as much as 1.23db. 

Similarly, the bandwidth shifts can produce effective 
cell length changes that result in actual cell area changes, 
consequently these shifts must also be accounted for in the 
evaluation of the cell area, A. This is automatically 
accounted for in the TAMU data reduction system because of the 
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way the area is calculated [2]. This is discussed further in 
Section S.O. 

4.2 Effects on Usage of Antenna Cain and Roll-Off Filter 
Funct ions 

The shift in filter center frequency must be accounted 
for in the system antenna gain, G r G t , end in the evaluation 
of the roll-off filter function, Z(f). The effect on roll-off 
was found to be of less consequence than in the antenna gain. 
Because of the reportedly irregular shape of the antenna gain 
pattern, small frequency shifts tend to produce substantual 

A_ 

changes in the resulting calculated value of o . The most 
dramatic effects were noted in the 13.3 GHz and 4.75 GHz 
antenna patterns. This effect is also explained in the last 
few viewgraphs in Appendix R. This effect has been subsequently 
diminished to some extent by using smoothed antenna gain 
tables [4]. 

4.3 Low Viewing Angle Effects 

The analysis of geometry effects on the solution of 
equation (1) reveal that there are limits on the doppler band- 
width at low viewing angles. As illustrated in Figure 6, 
whenever half of the actual bandwidth, BW^/2, is greater than 
the doppler center frequency, f dc- for a specified filter; 
e.g., -5 degrees viewing angle, the lower bound of the band- 
width will extend off the end of the positive spectral filter 
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look angle* ® f 
roll angle* <|> 



Isodopplers Corresponding 
to the Desired Bandwidth 
(BUD). If *F0 then f d < 
can occur. ' , 


Action To Be Taken: 

. If f dc < BW/2, set f dc • BW/2 

or 

• *4c < set * 


Figure 6. Bandwidth and frequency limits at low 
viewing angles. 


region. This has the effect of limiting the actual viewing 
angle to certain minimum values. It can he shown that these 
minimum viewing angle values are given by [2) 

k ■ Wsi* 

and 

L ■ ground cell length 
H ■ altitude of aircraft 
$ ■ drift angle of aircraft 

<j» • roll angle of aircraft 

Some values of the above expression are illustrated in Figure 7. 
Note that for increasing cell sizes the limit rapidly approaches 
-5° when only a small amount of aircraft roll is applied. 

Two approaches are used to handle the above limiting 
cases. At TAMU the value of the true viewing angle is limited 
so that the effective bandwidth never extends off the end of 
the PSD. At NASA/JSC the output data set is flagged with a 
caution note. Either technique is considered acceptable. 

5.0 GROUND RESOLUTION CELL SIZE 

Two different definitions of ground cell size or ground 
cell length were found to be in use by TAMU and NASA/JSC. 

First, TAMU software maintains a constant cell length. Figure 8, 


n 


Roll 




Figure 7. Viewing angle limits as a function of 
roll angle and drift angle. 
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Typical Cell Lengths -TAMU Processor* 



T|-T 0 (sec) 

L (meters) 

L(feet) 

L' (meters) 

L'(feet) 

13.3 GHz 

0.08 

25 

82.02 

18.68 

61.3 

1.6 GHz 

0.41 

50 

164.04 

18.42 

60.4 

0.4 GHz 

0.82 

75 

246.06 

11.84 

38.84 


♦Computed for 2048 samples/record and 150 Kts ground speed 



Bandwidth is adjusted to keep L constant for 
all and velocities. 


Figure 8- Illustration of the ground cell length as 
defined by TAMU. 
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as aircraft parameters vary. This cell length includes both 
the instantaneous viewing length and the "smear" distance 
moved by the aircraft during the integration time , t Q to tj. 

On the other hand, NASA/ JSC defines the cell length as the 
instantaneous viewing length, Figure 9. In the NASA/ JSC 
program the cell length, L' , is a program input, and bandwidth, 
BW., is adjusted to keep the input value of L' constant for 
all viewing angles. However, as the aircraft flight parameters 
change the total ground coverage per cell, L, computed using 
the NASA/ JSC processor will change slightly. This difference 
in the TAMU and NASA/JSC processors is not significant. But , 
the fact that NASA/ JSC defines the cell resolution length as 
the instantaneous resolution can be significant if the user 
does not understand what it means . It is important since the 
actual ground coverage is greater than the instantaneous 
resolution by the "smear" distance that results from the time 
integration process involved in filtering. 

6.0 COMPARISON OF GENERAL PROGRAM STRUCTURES 

The over all systems used to compute o° at NASA/ JSC and 
at TAMU are similar in the sense that both use the discrete 
Fourier transform to calculate the estimated PSD; however, they 
differ quite significantly in the way the data are handled 
before and after the PSD estimate is made. Also, there is a 
basic difference in the way the sign-sensing is accomplished. 
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The general approach used by both TAMU and NASA/ JSC is 
illustrated in Figure 10. Phase I consists of digitizing 
selected segments of analog data from the 14 -track tapes made 
on board the aircraft. Phase II uses the 9-track CCT containing 

a 

the digitized data to calculate and record o along with 
appropriate aircraft parameters and time. 

0.1 NASA/ESD General Program Structure 

The NASA/ESD general program structure is illustrated in 
Figure 11. Note that the system builds two complete and 
separate files, one containing ADAS/NERDAS data and one con- 
taining filtered data for each record from the CCT digital data 
tape. After the last record is processed the two files are 
merged into an output data set by averaging both aircraft and 
filter data over approximately 0.5 second intervals for 13.3 
GHz data and longer intervals for the lower frequencies. 

Next, these data are correlated to aircraft nadir time. The 
only time limitation noted in examination of this approach 
was that the two files built prior to merging dictate a limit 
on the length of the time segment that can be processed during 
any one execution phase. This limit is most severe for 
13.3 GHz data where only lines of 2S0 seconds or less can be 
processed. 

The sign sensing and PSD estimates are done by first 
calculating two discrete Fourier transforms, one for the sine 
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PHASE I 





Output 

Products 


Figure 10. Scatterometer data processing procedure. 
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Figure 11. General program structure of the 
NASA/JSC processing system. 



channel and one for the cosine channel of the input data set. 
The output of these two data sets are then used to get the PSD 
coefficients . for either the fore or aft data using the 


expressions 


m £(*S'Xt)\-(*c+nyJk' 

/t' = /c&+rc ) *V (KC-JS?'] A. *• 


where 


* fore PSD coefficient of i z frequency component 

■ aft coefficient of i frequency component 

t>h 

■ real Fourier coefficient of i frequency 
component, sine channel 

* imaginary Fourier coefficient of i th frequency 
component, sine channel 

“ imaginary Fourier coefficient of i*h frequency 
component, cosine channel 

t* h 

s real Fourier coefficient of 1 frequency 
component, cosine channel 

* scaling constant 


6.2 TAMU General Program Structure 

In contrast to the above, the TAMU system illustrated in 
Figure 12 processes each digital record completely, from 9- 
track input CCT to the output data set, in one pass. Using 
this technique there is no specific limits on how long each 
input data set must be. The processing begins with the first 


L 



r START J 


READ 

SET-UP 

DATA 


INITIALIZi 
SYSTEM 
CONSTANTS 
BUFFERS 


READ 

DIGITAL 
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COMPUTE 

f A* BW A’ 0 AL’ 
For each cel 


PERFORM CX 
-DFT, EVALUATE 
FILTER PUR 
& CALIB. PUR 


CALCULATE & 
BUFFER 
A/C DATA 
CORRELATED TO| 
imx 




START/STOP TIMES 
TAPE FILE # 

CONSTANTS, DUMP FUGS 




PRINT 
BUFFER 
LINE 



Figure 12. General program structure of the TAMU 
Processing System. 
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record and continues serially until the last record is done. 

The other important difference is in the way the discrete 
Fourier transform is used to do both PSD estimating and sign 
sensing in one pass. As illustrated in Figure 13, the NASA/ESD 
system requires two fast Fourier transforms to produce the aft 
data PSD, while with one complex fast Fourier transform the 
TAMU system produces both the fore and aft PSD. 

Actual data runs of the TAMU system have established that 
the average machine computation time per record is 0.102 
seconds. The average I/O time is 0.016 seconds per record for 
a total processing time of 0.118 seconds per record using the 
TAMU AMDAHL 470. Using the latest pricing equation at the TAMU 
Data Processing Center the cost per record comes to $0.01334. 

Note that it matters some if the data is 0.4 GHz. In 
this scatterometer system the sign sensing is not done, and 
thus the TAMU system, as with NASA/ESD, performs two fast 
Fourier transforms to get the PSD estimate. Run records for 
0.4 GHz data through the TAMU system indicate an average of 
0.1763 seconds per record with the I/O time the same; i. e., 
0.016 seconds per record, giving a total of 0.1923 seconds 
per record for 0.4 GHz processing. 

6.3 Comparison of Processing Results 

Comparisons were made of the output of the two scattero- 
meter systems. The sample data set chosen to use for the test 
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NASA/ JSC Double FFT 



Aft Spectrum 



Fore end 
Aft Spectrum 


Figure 13 Sign sensing techniques used by TAMU and 
NASA/ESD. 
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was Line S, Run 2, from Mission 347 made on 23 September 1976. 
Data acquired over Fields 17 and 18 were used as test data. 

In one test both TAMU and NASA/BSD processed the data using the 
same digital input data (Figure 14a) and b)). The digital data 
tape used for this test was produced by NASA/ESD. Some 
reformatting was done by both NASA/ESD and TAMU to make the 
input compatible with TAMU software. The results shown in 
Figure 14 indicates that both systems gave the same estimated 
o within +1.5 db. This test demonstrated that Phase II of 

m 

the processing systems were comparable. 

A second test was run at TAMU using the same data time 
segment, but processing was done using as input a digital tape 
made at TAMU. The results of this test, shown in Figure 15 a) 
and b) , indicate that both systems produce the same result from 
end-to-end. This test demonstrates that Phase I, the digitaza- 
tion process, is also comparable. 

Some problems in getting equivalent results were encoun- 
tered early in testing, but subsequently were resolved as 
differences in the system constants, antenna patterns and 
roll-off tables between the two systems. 

7.0 TAMU COMPLEX DIGITAL FOURIER TRANSFORM SYSTEM 

7.1 Advantages of New System Over the Hilbert Transform System 
The Complex Discrete Fourier Transform (CDFT) software 

system is designed to interface with the same digital computer 
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Mission 347 Line S Run 2 (Lawrence, Kansas) 
September 23, 1976 16:57:58,0 to 16:58:01.5 



-10 -20 -30 -40 


Look Angle (degrees) 


Figure 14. a) Comparison of TAMU and NASA/ESD processing 
software using as Input a common digital 
tape generated by NASA/ESD. 
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Figure 15. a) Mission 347 Field 17 processed by TAMU 
using as Input a digital tape qenerated 
by TAMU. 
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Figure 15. b) Mission 347 Field 18 processed by TAMU 
using as input a digital tape generated 
by TAMU. 
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respective filter, as opposed to an average center frequency 
used by the earlier HT system. Similarly, the value of the 
ground cell area. A, is calculated using the actual upper 
and lower doppler contours that bound the respective cell being 
viewed at that instant in time, rather than average frequencies 
for each cell over the flight line. 

The ratio of the power received to the power transmitted 
is taken from the sum of the power spectral density (PSD) 
lines of the CDFT within the doppler frequency band defining 
the ground cell can be continuously updated as a function of 
aircraft parameters, the DFT filtering approach provides a 
technique to adaptively filter the radar output to maintain 
cell size and incident angle. 

In the HT system the digital filter is used to acquire 
power within each band, including the calibration power, P t< 
Each of the digital filters is pre-designed for an average 
center frequency over the entire flight line. This does not 
provide constant resolution or incident angle as aircraft 
flight parameters vary since aircraft dynamics cause the 
center frequency corresponding to each viewing angle to vary 
over the length of the flight line. To calculate a new center 
frequency and re-design each filter for each ground cell 
would be absolutely prohibitive in machine costs. However, 
the nature of the CDFT makes it quite easy to optimize band 
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compatible tape (CCT) format used by the Hilbert transform (HT) 

sysyem [1] previously developed at TAMU to process scattero- 
meter data. The system performs repeated scattering coeffi- 
cient estimates for eight aft viewing angles. The output data 

A- 

set has the scattering coefficients, a , aligned with respect 
to the aircraft nadir times. As before, the scattering coef- 
ficient is calculated using the expression 


A 


Q/rrf Jk 

y? 4 JCs 


where X * wavelength 
R = range 
A = area 

= cable loss 
K g = system constant 
Z(f d ) = roll-off 
G r G t = antenna gain 
P r * power received 
P = power transmitted 


( 16 ) 


The important difference in the solution techniques used by 
the CDFT is in the evaluation of the four variables R, A, P , 
and P t « 

First, the range, R, is calculated to the center of the 
actual ground cell based upon actual center frequency of the 
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and center frequency for each cell evaluated by the stem. 

7.2 CDFT Algorithm Summary 

For each set of digital data taken from the CCT, the CDFT 
system must evaluate -sight sets of range, area, roll-off, gain, 
and power-ratios. This is accomplished in the order illustra- 
ted in Figure 16. The symbols used are defined in Table 1. 

To evaluate 6^, the doppler angles, the viewing point coordi- 
nates of each viewing angle are first calculated using the 
geometric relations in Figure 17. 


y' . //tW V 

(17) 


(18) 

X - X 'ces $ + y's/»j6 

(19) 

y r X'^^0 - y' Cos X 

(20) 

The doppler angle, 0^, is then given by 

-,r X* 7 ^ 

3/ 

(21) 

The doppler cent* . r frequency is then given by 

£ r ( 

(22) 
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Figure 16. Computational algorithm for computing 


sigma zero. 
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TABLE 1 


Symbol Cross Reference/Definition 

0* Doppler Angle - Viewing angle referenced to the plane of 
the ground track and aircraft gravity vector. 

f. Doppler frequency associated with a particular viewing 
angle. 

Af d Bandwidth desired or calculated. 

NF- Number of filter elements (or discrete spectral lines) 
needed to make Af^. 

BW t True bandwidth - NF^ times the spectral line width. 

N c Center frequency index value. 

f. True doppler (center) frequency associated with the 
true bandwidth, BW t< 

*fdl Left-most doppler filter index pointer 

Ifdr Right-most doppler filter index pointer 

e dt True doppler angle associated with true doppler center 
frequency, f d(; . 

R t True Range - range from antenna to the true center of 
the ground cell defined by f dc and 0 dt . 

e A Viewing angle referenced to the aircraft coordinate 
system, but without pitch correction, used for gain 
look-up. 

0,. True viewing angle, corrected for shift in doppler 
center frequency due to discrete PSD lines. 
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to - roll angle 
<p - drift angle 

Figure 17. Definition of processing coordinate system. 
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where 


v ■ aircraft ground speed 
X * scatterometer wavelength 


Next, the desired bandwidth, Af^, is calculated using 

A 6 = Xk.to* S 6<* 


( 23 ) 


where 


Xk. - Cx) H 


( 24 ) 


and L » instantaneous ground cell length 

H ■ aircraft altitude 

The number of filter elements required in the PSD set to repre- 
sent the desired bandwidth is calculated by 


A/?i 



( 25 ) 


where N ■ number of lines in the PSD. 

f * sample frequency of digitizer when making 
s the CCT. 

The notation [ ]j means integer value of the floating point 
expression enclosed. 

Note that the true bandwidth is the sum of the filter 
elements needed, NF^, thus 



gives the actual bandwidth over which the power will be sumed. 
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The center frequency index for the PSD buffer is given by 


A/c -- [&(%)+* €* (27) 

and the true center frequency is then 

& - A/c. (ty*)) C2») 

whenever the value of NF^ is odd. When the value of NF^ is 
even N c is given by 

t29 ’ 

and the true center frequency is 

^ =. (a/ c (30) 

Similarly, the value of the PSD left and right power sum- 
mation pointers depend on whether NF^^ is odd or even. For the 
case where NF^ is odd the left pointer is given by 

jr& t = -0/*- (31) 

when NF^ is even, the left pointer is given by 

* A/c, A/ - A/fi/z. (32) 
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In either case the right pointer is given by 


Next, since the actual doppler center may have been 
shifted by as much as one-half a spectral line in the PSD, the 
true doppler angle, and thus the true viewing angle, may be 
shifted accordingly. The true doppler angle is given by 


* STAJ 



( 34 ) 


The true (or corrected) viewing angle is now evaluated by 

[ # 


7 


( 35 ) 


The true range, R t , to the cell may be evaluated by 

ZJ - ///<**!&£ (56) 

and the antenna angle, d^, is given by 

0 A , 'A t ( 37 ) 

The value of antenna gain is evaluated using the actual 
viewing angle, 0^. The value of system roll-off is evaluated 
using the true doppler center frequency, f d( .. 
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The value of cell area, A, is calculated using the doppler 
contour technique described in [4], [5]; i.e., 


A* 


to 'Cy± - y*) 


<fos 


where 


1 ^' - 2. A/m. 

(fos y* 


antenna beamwidth 


The value of y ^ and y 2 are given by 


Yi - [ *“*** r + (*»**(**->)+ ki<*» V-/) K J 


where k. 


( rr ) 2 


The value of f^ is interpreted to be f^ * lower doppler fre- 
quency limit of the actual bandwidth and f 2 ■ upper doppler 
limiv of the bandwidth. 

After completing the above calculations for each of eight 
viewing angles, the PSD buffer is calculated using the CDFT. 

The input to the transform is programmed so that the real part 
is always the true sine signal and the imaginary part is always 
the true cosine signal from the radar analog data set. The 
output of the CDFT, which convolves on the basis of e is 
the set of Fourier coefficients arranged such that the aft 
data is associated with positive frequencies and the fore data 
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is associated with the negative frequencies, Figure 18. The 
PSD buffer is then filled with the square of complex absolute 
values of the aft data coefficients. 

Using the index pointers already computed, the P f and P ( 
values are found by summing over the appropriate sets of PSD 
lines. Now one only needs to evaluate equation (1), above, 
for each of the eight viewing angles. This will complete the 
calculation cycle for one instant in time, and one CCT record. 

The TAMU CDFT system is designed to process data from any 
one of four scatterometer frequencies; i. e., 0.4, 1.6, 4.75, 
or 13.3 GHz. The system has been tested extensively using 
1.6 and 13.3 GHz data inputs. It has had some testing with 
both 0.4 and 4.75 GHz data, however, measurements at these 
frequencies were limited during the course of this evaluation. 

Appendix C is a FORTRAN listing of the processing system 
and is complete with the exception of the antenna gain tables. 

Appendix D contains references [3J and [4] which might 
otherwise be difficult to obtain. 
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Figure 18. Spectral arrangement of the complex 
digital Fourier transform. 
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APPENDIX A 


FLOW CHARTS FOR HILBERT 
TRANSFORM PROCESSING SYSTEM 
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APPENDIX B 


VIEWGRAPHS OF SEPTEMBER 1977 


(J 


SCATTEROMETER DATA PROCESSING 

Summary 

Discrepancy: Misinterpretation of aircraft flight logs/tape I.D. 

1. IS dB gain assumed in one channel of the 
1.6 GHz system. 

2. Misinterpretation as to what roll-off filter 

WAS USED ON SOME DATA RUNS. 

3. Terminology denoting in-phase and quadrature 
CHANNELS OF 1.6 GHz AND 13.3 GHz SYSTEM 
CONFUSING AND MISUNDERTSTOOD. 

Action: 1. Verbal clarification of meaning of flight log 

INFORMATION (CORRECTED 1. AND 2.) 

2. Spectrum analysis to determine channel containing 

CALIBRATION SIGNAL AND TIME HISTORY PLOTS OF SEVERAL 
ANGLES OF NON-SHIFTED PROCESSED DATA TO DETERMINE 
WHICH OF THE QUADRATURE CHANNELS LEAD AND WHICH 
LAGS IN-PHASE. 

Result: The calibration signal was found to be on opposite 

QUADRATURE CHANNELS FOR THE 1.6 GHz AND 13.3 GHz RADARS. 




NASA Tape Duplication Facility (Bldg. 12) 



iscriminato: 














SCATTEROMETER DATA PROCESSING (CONTINUED) 

Summary (continued) 

Discrepancy: Inconsistancy of time reference 

* 

j . 1. ADAS INFORMATION COULD NOT BE DISCRIMINATED 

• Poor tape reproduce capability (6FE recorder) 

• ADAS DISCRIMATOR ADJUSTMENT 

* 

2. IRIG TIME USED TO PROCESS DATA 

• Did NOT PROVIDE ONE TO ONE TIME CORRELATION 
BETWEEN PADAR DATA AND PHOTOGRAPH 

• Did NOT PROVIDE TIME CONS I STANCE BETWEEN RADARS 
EVEN WITHIN SAME RUN 

j Action: 1. Adjustments made to ADAS discrimator 

\ - 2. Better quality tape recorder obtained 

» 

j 3, ADAS TIME reference used for all processing 

f 

| . This provides time consistency between all 

t 

1 RADARS/ DATA RUNS/ AND PHOTOGRAPHY 


j 

1 


7 * 



SCATTEROMETER DATA PROCESSING (CINTINUED) 


Summary (continued) 

Discrepancy: Misalignment of sigma for each ground cell. 

1. Average flight parameters used to compute 

AN AVERAGE TIME SHIFT FOR EACH ANGLE DURING 
EACH RUN 

2. The average velocities that were used to compute 

SHIFTS INCORRECT. 

• AVERAGES WERE COMPUTED FOR TIMES ON 
FLIGHT LOGS 

• TIMES ON FLIGHT LOG INCLUDED SET-UP FOR 
DATA RUN AND USUALLY CONTAINED LARGE 
VELOCITY VARIATIONS 

• VELOCITY ON FLIGHT LOG NOT ALWAYS 
REPRESENTATIVE OF VELOCITY DURING DATA 
RUN 

Action: Program modified so that flight parameters are updated 

CONTINUOUSLY FROM ADAS DURING SlGMA COMPUTATIONS. TlME 
SHIFTS, AREA CALCULATIONS, ETC. ARE NOW AS ACCURATE AS 
POSSIBLE. 




SCATTEROMETER DATA PROCESSING (CONTINUED) 



r 

v 
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i 

•i 


Summary (continued) 

Discrepancy: Calibration Constant Used In Processing Data 

1. Compensation of calibration constant for 

THE MEASUREMENT TECHNIQUE IMPROPER 

2. Compensation of the calibration signal for 

SIGN SENSING IMPROPER (IN ANALOG PROCESSING) 

| 

5 

i 

Action: Program corrected after agreement was reached on I 

PROPER COMPENSATION j 

\ 

i 

1 

(Note that this discrepancy only introduced a constant offset 

IN THE PROCESSED DATA,) j 

j 

-i 

I 

4 

! 
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SCATTEROMETER DATA PROCESSING (CONTINUED) 


1 
■ > 

It 


Summary (continued) 


Discrepancy: Scatterometer hardware inconsistancies 


r 

! 




1 . 


1.6 W AND 1.6 VH DATA NO GOOD BEFORE 
MISSION 362 

• SWITCH IMPEDANCE MISMATCH 


2. 400 MHz HORIZONTALLY POLARIZED ANTENNA HAD 

A LOOSE CONNECTOR DURING JSME FLIGHTS 


Action: Corrected by E&D Personnel 


t 

i 

! 



SCATTEROMETER DATA PROCESSING (CONTINUED) 

SllffllA&X (CONTINUED) 


Discrepancy: Inability to reproduce raw data prom analog 

TAPES WITHOUT DROPOUTS, CLIPPING, OR ADDITIVE NOISE 

1. Analog GFE tape recorder available to TAMU 

WAS AN OLD FLIGHT RECORDER THAT HAD VERY POOR 
REPRODUCE CAPABILITY 

2. Analog data tapes are duplicates of poor 

QUALITY WITH CARRIER LEVELS EXTREMELY LOW 

3. Duplicate flight tapes given to TAMU for pro- 
cessing ARE OF POOR MECHANICAL QUALITY (REUSED 
TAPES) 

NOTE: Little quality control is apparent in 

DUPLICATION PROCESS 

A. FM MODULATION IN RECORDING RADAR DATA ON 
AIRCRAFT EXCEEDS THE A0% DEVIATION LIMIT 
OF THE RECORDERS MUCH OF THE TIME 

5. CONSISTANCY OF AIRCRAFT RECORDER SET-UP IS 
POOR, ADDING TO DIFFICULTIES IN 3. 


SCATTEROMETER DATA PROCESSING (CONTINUED) 

Summary (continued) 

6. Recorder calibration levels put on data tape 

ARE OF VERY POOR QUALITY 


Action: 1. A better quality GFE laboratory tape recorder has 

BEEN PROVIDED TO TAMU AND INTEGRATED INTO THE 
DATA PROCESSING SYSTEM. 

2. TAMU DIGITIZER MODIFIED TO BE COMPATIBLE WITH 
INPUTS OF ♦ 5 VOLTS IN ORDER TO HANDLE THE OVER- 
MODULATED RADAR DATA 

• Tests were performed to insure that the 

REDUCTION IN DIGITIZATION RESOLUTION DID 
NOT AFFECT THE FINAL SIGMA 

3. Quality control of data recording, tape reproduce 

AND TAPES MUST BE DONE BY NASA. 




IIICRQHAVE PROGRAM 




SCATTEROMETER DATA PROCESSING (CONTINUED) 

Summary (continued) 

Discrepancy: Uncertainty in antenna patterns 

1. 13.3 GHz data have consistent "humps" 
at 30° and 50° 

2. 13.3 GHz pattern (measured in 1976 by NASA) 
HAS STEEP GRADIENTS AT 30° AND 45* - 50° 

3. 1.6 GHz DATA HAVE CONSISTENT "HUMP" AT 
20 ° 


Action: 1. Attempt to approximate more realistic pattern 

2. Process data at 5°, 10% 15°, 20°, 25°, 35°, 
45% and 55°. 
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SCATTEROMETER DATA PROCESSING (CONTINUED) 

Status 

1. Currently processing portions of FCF over Death Valley to 

FURTHER INVESTIGATE ANTENNA PATTERN DISCREPANCIES AND AS 
FINAL CHECK ON SOFTWARE SYSTEM 

2. A LIST OF SPECIFIC ADAS TIMES ARE BEING PROVIDED TO TAMU 
BY NASA FOR PROCESSING 

3. Processing will commence after resolution of antenna 

PATTERN UNCERTAINTIES 

4. Character of 400 MHz data unexplained and must be taken 

AT FACE VALUE UNTIL MORE DATA ARE PROCESSED 
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LISTING FOR PROGRAM SCATTER.CZTA78 
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>** TEST PROGRAM FOR SIMULATING CHIRP-Z ALGORITHM *=< 

OEVELOPFO AT TEXAS A&M UNIVERSITY 
REMOTE SENSING CENTER 

SSPT78 • MAY79 

ALGORITHM PERFORMS COMPLEX OFT <X*JY) ON t V0 INPUT CHANNELS 


SYSTEM PROCFSSES SCATTEROMETER DATA FOR 13# 3* 1.6, 4.75* OR .4GHZ 
»»> VERSION OATEO MAY 1979 ««« 

***** input parameters: 

CARD *t POINTS* COL 1*4* 14 (1024 OR 2048) 

SVSK* COL 6-10* FS.O 
CARL* COL 11-15, PS*0 

DELTs COL 16-20. F5.0(NBR SECS TO PROCESS) 

TSC I ) * COL 21-22, F2. 0( START OF 

TSC 2 ) * COL 23-24* F2.0(TIME SEGMENT * 

TS(7)* COL 25-27, F3. I HHMMSSS ) 


CARD*2 85 CHARACTERS OF COMMENTS 

CARD*3 INCR* COL i-5,I5(EPT LENGTH) 

ITYPE * COL 6-10. IS (1*13.3, 2*1.6. 3*4.75, 4*0.4 GHZ ) 

DT s COL 11-1*. FS. 3, ( AVG RCD TIME INTERVAL) 
l BUG * COL 16-20 .15, (DATA DUMP SWITCH) 

KPOLZ* COL 21-2E. IS. (SET IF V-POLARIZED C-BAND DAT*) 
I PLOT s COL 26— TO, IS. (SET IF SPECTRUM PLOT OF 1ST RCD) 
CELL* COL 31-75. FS.0( RESOLUTION, INSTANTANEOUS) 


01 MENS t ON ZV<0) ,A(S).90W<9),!F0L<8).IFDR(8).GAINSQf 8),BEAM(S) 
DIMENSION MATRtXC19.100).R<8) .TSC3).FOOP(8).THETA(S).DFDOP(A) 
DIMENSION SIGMA(8) , I BPTR< 8 ) ,FREQP< 1024 ) • AFT( I 024 ) , AFTC( 10 24) 
DIMENSION TRIX< 16) , PSPFCf 8 ) , ANGTf 8) »PROPT( 8) , I THETA ( B) 

I N T EGER*2 TtME( 3» .ADAS (5). SINE! 2048) .COSINE ( 2040 ) , IC0M4 40 ) 
INTEGER POINTS 
REAL NTS 

REAL L AMD A »KONE » K T WO 

COMPLEX FRS( 204 8 ) » CM®L X , FRC( 2048 ) 

COMMON /FILT/TTVPF. KPOLZ 

COMMON /I TAPE/TJ ME. ADAS . S IN*. COSINE 


nn no on on on 


COMMON/KARE A/KONE* KTMO 
C3MM0N/RAT/F INCR *PCAL *FSAMR • f NCR 
OATA fORP!/I2»S6«3V 
DATA MATRIX/I 8004-999/ 

C AFT ANGLE SET IN DEGREES* UR TO EIGHT VALUES 
DATA THETA/S#*I0#*1S#*20#*2«#*3S# .40*. 45./ 

OATA I THETA/5* 10* IS* 20* 29*35 *40*45/ 

GFT SYSTEM RUN RARAMTERS AKO COMMENTS 
REA0I5 *1 0 1 R01NTS*SYSK*CPBL*0ELT*TS 
10 FORMAT! I 4*1 X*3FS#9*2F2*0*F3#1 I 

COMMENTS INCLUOE MSN.LINE *RUN*D4Tf *8AND IOENT.ET# 
READ! 5*1 1 IICOM 
It FORMAT! 40 A2 I 

REAO OFT LENGTH* TYRE SCAT* AND AVG RECORO TIME DEL . A 

SET tSUG SVIYCH IF OATA DUMP OFSIREO ON UNtT 9 
READ! 5*121 INCR* ITYPE*0T * IOUG*KROLZ* tRL0 T *CELL 
12 FORMAT!21S*F5#3*3IS*FS#OI 
GO TO ! 1 *2*3*41 • I TYPE 

1 3* 3GHZ SCATTEROMETER 

1 LAMDA*0#022S4 
FS4MP=2500 r * 

ECAL*1 2090* 

FNOI 7=8000* 

BWI0*6#3 
GO TO 13 

C 1#6GHZ SCATTEROMETER 

2 LAM0A*0*1GT3T 
FSAMP*5900* 

FCAL«1909# 

FN0IZ*1S*0* 

BW l 0*7# 94296 
GO TO 13 

C 4#75GHZ SCATTEROMETER 

3 LAMDA=0* 06311 
PSAMR*10000# 

FCAL*2930# 

FCALV=3380# 

IF«POLZ#NE*0l ECAL*FCALV 
FNOI Z* 31 00* 

BWI 0*2*5 
GO TO 13 

C 0#4GHZ SCATTEROMETER 

4 LAMOA« 0*74948 
FSAMP*2500# 

FCAL*1000# 
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FCPT*500* 

FNOIZ* 200* 

B»IO*14«0 

C 

13 FPNTS*FLOAT( POINTS) 

TM5T8»TANf TAI 8»/ST« 1 » 

ML*100 
I CLR*0 
PT4»1# 

PC4*1 • 

C SET UP DEFAULT VALUES FOP CELL LENGTH 
IF(CELL«GT«0*> GO TO IS 
GO TO <101 *102* 103* 104 ) • ITVPE 

101 CELL*25« 

GO TO 15 

102 CELL=50* 

GO TO 15 

103 CELL»40* 

GO TO 15 

104 CELL*75« 

15 CONTINUE 

C TO CORRECT FOP MEASUREMENT OF CALPEF 
IF! I TYPE* EQ# 2) SVSK*SVSK4'5#0103*2# 

VK0B«30«*AL0GI0f FOPPI >-20.* ALOG1 KLAMDA) -SYSK+CABL 
C POWER SPECTRUM LENGTH * « LENGTH Q* FFT)/2 # 

FlNCR*FLOA T ! INCP* 

INCR2*INCR/2 
DELFPQxFSAMP/lFINCPI 

INITIALIZE OUTPUT LOAD POINTER 6 LINE COUNTER 
MXRCTRsML-10 
NOUTsO 

C NTS IS THE INV*»SE OF SPECTRAL LINE W JDTH!L INES/HZ 1 
NTS»FI NCR /f SAMP 

C CALCULATE POINTERS to CALIBRATION VALUES IN AFT DATA BUFFER 
MC*IF!X<FCAL*NTSO*5)M 
MC4*IFIX<FCPT*NTS*0*S)41 

C CALCULATE POINTFRS FOP NOISE BAND MEASUREMENTS 
IFNOIZ*IFIX<FNOIZ*NTS*0.5> 

NOIZL*lENOIZ-2 

NOIZRaIFNOTZT?, 

TSTART *3600* * r S! 1 ) A69»*TS ! 2 )*TS<3 > 

C CALCULATE NBR OF PCOS TO PROCESS * ! TIMF INCREMENT > /! RCD LENGTH) 
TS T OP*OELT/OT 
I END*I FIX! TSTOPO#B) 

NR t TE (6*60) ICOM*POtNTS*LAMDA*FSAMP*FCAL •SVSK»CABL*DELT »TS*IEND 
69 FORMAT! I HI »5X .4 9A2 • // • 1 SX • • POINTS/RCD** • IS* 5X • • PAVE LENGTHS • ,F7, 
SSX* ' SAMPLE FREO.*» »FT.0»5X* »CALIB SIGNAL FREQ**.F6.0./. 
$15X.»PEFFR CAL IBS* »F6# 1 * • OB* .SX .• CABLE LOSS®* *F6. 1 .SXt / % 

$ 1 SX* *NBR SEC PROCESSED** »F5« 1 »SX» •START TIMFS • • 2F 3. 0 .F4. ! , 
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$SX.*NBR RCDSl * . I6./I 
WRITE! 6.16)NTS«MC. IFNOI Z 

16 FORMAT! 1H .1 AX, 'SAMPLE INTERVALS • .F9. 5.5X. *CALIB PWR PTR:*.tft. 
S5X.*N01SE PWR PTR:*.I6«//I 

20 CALL INPUTIPOINTS.TI 

C DIGITIZER CHANNEL 1* StN. CHANNEL 2* COSINE. IN CORE 

C IF 1 3. 3GHZ TAPE CHlt *> SINE. ▼APE CHI 3 *> COSINE. ( CHI 1 *TS I NE. CHI 3a 

C IF I.6GHZ TAPE CM 3 «> SINE. TAP* CH 1 *> COSINE. !CH3*TC0S.CH1*TSI 

C IF 1 • 6GHZ TAPE CH T a> SINE. TAPE CH 5 »> C0SINE1CR0SS POLZI 

C IF 4.75GHZ!LtKE) TAPE CHlaTSINE. CH3«TCOSINE 
C IF 4.7SGHZ! CROSSI TAPE CH5*TSINE. CHOaTCOSINF 
C 

IF«TSTART,GT.T| go to 80 
I Ft MOOT. EQ.O) WRITE! ft. 81 I TIME. T, ADAS 

21 FORMAT 1 1 HO *5X. • F IRST PECORO TIME • *2( 12. • ! • » .13. 

*10X.»TIME IN SECSS* ,F**« l ,/ , t OX. • ALT: • • 14 . 

S5X.* VELOCITY:* .I4.SX. • pitch: • .I4.SX.*ROLL: • .14, 

SSX**DRIFT:« .14./) 

HaADASU I 

H*H*O.3O40 

IFIH.LT. 350. 1 H*450* 

VEL*ADAS!8) 

VEL*VEL*0«514 

|F( VFL.LT.60. * VEL*TT. | 

C AOA S/NERO AS FORMAT ON TAMO OtGITAL TAPES HAS PITCH FIRST 
PITCHaAOAS<3) 

PITCHaPI TCH/5T3. 

ROLL*ADASt 4 1 
POLL “ROLL/ST 1 
ORtFT*AOAS<S) 

DR IFT =DR IFT/5T3, 

FILTER THE AIRCRAFT DATA FOP OLO OtGITAL TAPFS 

LIMIT * 5 DEG for PITCH. 10 OEG FOR ROLL. AND 15 OEG FOR DPIFT 


API TCH*ABS< PITCH ) 

A ROLL* ABS! ROLL I 
ADRIFT«ABS(OR TFT ) 

IF(AP!TCH.GT.0.0S?86t PITCH»0.08T£6*< PITCH/API TCH I 
IFf AROLL.GT.0.1T4S8) RCLL*0. IT458*«R0LL/A90LLI 
IF! ADRIFT.GT.O.261T01 DPI FTsQ. 2ft I 78 *( D»! FT/ADR IFT ) 


CALCULATE NBR OF CELLS IN VIE* 
XCELLSsH*THET9/C VEL*OT) 

INITIALIZE PRINTER -JINTFR 81 AS 

IFtNOUT.EO.OI LPRNT*IFIXf XCELLS*! .3) 

COMPOTE SIZE OF INSTANTANEOUS RESOLUTION CELL 
SL I P*V8L*FPN T S/FSA MP 
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SCELL*CELL-SLIP 

COMPUTE BUFFER LENGTH NEEDED 
MLN*IFIX< XCELLS*0*5I 
IFf MLN«GT*MLI GO TO *00 

calculate theta function data values 

VEL20L*2«EVEL/LAMDA 

XK*SCELLEVEL20L/H 

IF<N0UT*EQ*0I WRITE! 6*221 H ,VEL. PITCH, ROLL .DRIFT ,SLI P. XCELLS . 
* SCELL • BW t O • CELL » I NCR • I TYPE • VEL 20L 
22 FORMATIIH .SX.'RUN PARAMETERS* ,/,6X, *H5 *.FS*1, 

*5X,*VFL: • ,F6#1 ,6X»» PITCH S' ,F6* I »SX» 'ROLL S • #F6* I • 

S5X. 'DRIFTS • ,F9. t,/,6X, 'SLIPS ' »F6. t ,3X, 'XCELLSS' ,F6. I , 
S3X.'SCELLS' ,F6*1,SX.'BEAMWS' »F6*l *2X*'CELLS' •F6*I*5X« 

S'lNCRS I5*/.6 X,'ITVPE:' ,IA*«x,'2V/LAMOA:' *FR # 2,/I 

DO 34 I*t.S 

FOR EACH AFT ANGLE IN THF SFT 
ANGL*-TMETA< l I /ST* 3 

GUARD AGAINST DOPPLERS OF NEGATIVE FREQUENCIES OFF FILTER BANO 
AANGL*ABS< ANGLI 

IF ( AROLL*GE* AANGL I ANGL**AROLL 

COMPUTE VIEWING POINT COORDINATES ON THE GROUND 
VP I *««H*T ANC ROLL I 

XPI*SQRT| |H*TAN( ANGLI >**2*VPI««2) 

TTPa-XPI *COS< ROLL I /H 

XTs-H/COSIROLL »*ttp 

VT*H*TAN< ROLL ) 

XI = XT*COS< DRIFT I ♦V T PS I NI DRIFT | 

VlsXT*SIN<DRlFT»-YT*COS<DRtFT| 

FORCE THE DOPPLFR BAND TO BE ON FILTER BANK 
IF( XI *LT , SCELL/2* I XI*SCELL/2* 

COMPUTE CORRECTEO DOPPLER ANGLES 
XN*XI**2 
XO*H**2«VI**2 
ANGLO* AT AN! SORT! XN/XOI I 

COMMUTE DOPPLER CENTER FRFOUENCV «»■ ■ »■«■■■» ■ 

FO0P( I l=VEL20L*SIN! ANGLO I 

IF! IT YPE*EO* ♦ I FOOPI I l*SOO*-FOOP< II 

COMPUTE DOPPLER BANOWIDTM mm mm mm 
DFDOP! I l*XK*COS( ANGLO I E* 3 


COMPUTE NBP OF FILTER ELEMENTS NEEDED 
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NPt»lFtr.«DFDOP<t»*NTS*O.S> i 

1P<NFULT«I» NFI*» 

CALCULATE TRUE BANDWIDTH 5 

DFOOPI I »*FLOATI NPt »*OFLFRQ | 

IFtHOD(NFt»2t.EQ*0) GO TO 25 

ODD NBR OF FILTERS REOUTRED j 

CENTER FREQ INDEX ISt \ 

NC*lFlXfFDOPtI »«NTS+0*S> j 

TRUE DOPPLER CENTER FREQ 
FDOPI I )«NC*DELFRO 

LEFT INDEX POINTER 

IFDLI I l»NC-«NFI-| 1/2 

RIGHT INDEX POINTER 
GO TO 27 

25 CONTINUE 

EVEN NBR OF FILTERS REOUIREO 

CENTER FREQ INDEX 

NC*!FIX(FDOPm*NTS> 

TRUE DOPPLER CENTER FREQ 

PDOPII »*fFLOATINC»*0.5l*OELFPQ 

LEFT INDEX POINTFR 

IFDLCI »*NC+l-NFI/2 
2T CONTINUE 

IFf IFDL( I I.LE.O) IFOL«I)*1 
RIGHT INOEX POINTER 

IPDRI II*IFOLIl l*NFI-l 
IF| ITYPE.EO.AI FDOP< T >*SOO.-FDOP(I » 

COMPUTE TRUE DOPPLER ANGLE AND VIEWING ANGLE 
ANGTDsARS INI FDOP< I ) /VEL2CL ) 

XTSQ*XD*TAN< ANG t D)** 2 
VTSQsT 1**2 

ANGTL* AT AN< SORT ( X^SO+VTSO ) /H » 

C CHECK FOR NEGATIVE OOPPL»»S 
AROLLsABSCPOLL i 

ifiangtl.lt, arolli angtl*aroll 

C SET UP FOR AFT VIEWING ANGLES ONLY 
ANGTf I )»ANGTL*57.S 
C COMPUTE TRUE RANGE 

R( I ) *H/COS( ANG T L I 
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C CALCULATF AREA VARIABLES KONE A KTWO 
ICON*! 

KONE*A*/(LAMOA*(FDOP< I 1-0F00PC I »/2* ) )**2 
C CHECK POR ZERO S NF6A T t V? OOPPLERS 

IFf t OFOOPf t 1/2* I •GE.FOOPC 1 1 » I CON*- 1 
KTWO*A* /<LAMOA*I FOOP< I )tOFDOPI I )/2. 1 1**2 

•*> THIS SYSTEM FOR AFT ANCLES ONLY ««<****** 

COMPUTE CORRECTED VIEWING ANCLE THRU ANTENNA PATTERN COORD SYSTEM 
ANCLPI*ATAN<-COS< ROLL )*SORTITAN< ANGTL )**2-TAN<R0LL)F*2l ) 

• FIND VALUE IN TABLES OF GAINSO 
GAINSOI I )*GAI N( ANGLPT .PITCH) 

FIND VALUE OF WATER FILTER ROLL-OFF 
ZW< t )«GAMMA(FOnP< I > ) 

CALCULATE CELL AREA BETWEEN ISOOOP LINES USING TPUE LOOK ANGLE 

BEAMII)*BWro/ST,<5 

CALL CAREAIA< I I.BEAMf D.H.ANGTL.VEL.ROLL.ORIFT, ICON ) 

30 CONTINUE 

SET UP TO 00 r HE FFT WITH • INCR* PCtNTS 
NU*IFIX< ALOGIFfNCR >/AL0G<2« )TO*S) 

CONVERT DATA TO FLOATING PCINT 6 PLACE IN FFT BUFFER 
GO TO 131 .32,31 • 32 > . I TYPE 
32 CALL SCALXf POINTS. SINF .COStNE.FRS) 

IF«ITVPE*EO*A) CALL SCALXC POINTS. COSINE. SINE. FRC ) 

GO TO 331 

IF 13.3GHZ OR A.TS DATA < TAMU INPUT TAPES) PUT CHAN2 DIG DATA FIOST 

31 CALL SCALXIPOTNTS. COSINE .SINE. EPS) 

331 CONTINUE 

IL*POt NTS/I NCR 
DO 35 1*1 .IL 
IS*C 1-1 )*INCPT1 

BUILD SPECTRUM FOR <CHl)+j(CH2) 

CALL DFT < FRS( IS) .NU. I NCR I 
BUtLO CALIB CHANNEL SPEC FOR 0*AGHZ DATA 

IF< ITyPP.EO.AI CALL DFTIFRCf IS) .NU.INCR) 

3S CONTINUE 

SUM THE IL SETS FOR AFT DATA 
DO 36 1*1 . INCR2 
FREOPI I I *FLOAT< 1-1 )*OPLFFQ 
AFT ( I ) *0*0 
AFTCC I )*0.0 
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BUILO AFT DATA POWER BUMM OF • INCA/2* SPECTRAL LINES 
00 37 K«I*IL 
IS«!K-1)*INCR*I 
AFT! I >*AFT! I l+CABS! FRS ( tSI I 
IF|ITVPS«NE*4» GO TO 77 
AFTCIt )*AFTC!tl*CABS!FRC!tSlt 

37 CONTINUE 
36 CONTINUE 

IP(N0UT*GT»0t GO TO 74 

CHECK PLOT FLAG* BTPASS IF ZERO 
!F(IPL0T*E0«0) GO TO 74 
WAITE! 6*33 I I CON 

33 FORMAT! I HI *//*SX** AFT DATA SPECTRUM FOR PFC0»O NBR t • *SX* 40A?, /» 
CALL PL0T|FRE0P.AFT.tNC«2l 

tF!ITVPE*NF*4| GO TO 34 
WRITE!6* 1 33) JCOM 

133 FORMATOHl .//•S***AFT DATA SPECTRUM FOR 0*4GHZ CH2* BCD IS 
S5X*40A2*/» 

CALL PL0T!FRE0P.AFTC.INCR8» 

34 CONTINUE 

GET VALUES OF NOISE POWER 6 CALIBRATION POWER LFVELS 
CALL PWR!NOIZL*NOtZR»AFT*FNOm 
FNOIZ*FNOIZ/5 
SUN CALIBRATION POWER 
LPP*MC«6 
LRP*MC+6 

CALL PWR!LPP*LRP*AFT*PSU»T» 

IF! !TVPE*NE*4) GO TO 39 

SUN CAL IR CHANNEL POWER FOR 9*4GHZ CM2! AT 1 KH7I 
CALL PWR!LPP*LRP*AFTC*PC4» 

LPP4«MC4«6 

LRP4*MC4+6 

SUN TRANSMIT POWER FQR 0*AG HZ CHZIAT 0*5 KHZ) 

CALL PWR!LPPA,LRPA*AFTC,PT4I 

39 CONTINUE 

!F(NOUT*GT*OI GO TO 78 
WRITE! 6 *1 38 ) FNO I Z • PSUBT . PC4 . PT4 
136 F0RNATI1H **FNOtZ: • * fo c s*sx*«psubt:« *F9*S* 

»5X.»PCAL-4t* *FQ,S.SX.*PT«i4: •*F9*9» 

38 CONTINUE 

BUILD SIGMA-ZERO S^T 
DO 40 1*1 *8 

FIND VALUE OF RECEIVED POWER 




CALL PBRCIPOLC M.IPORC I> .AFT.PSUBRl 

SAVE DC BUG DUMP VALUES 
PSPECC I I *PSU8R 

PROPT < I >«10»*ALOG19C C PSUBR/PSUBT)»CPC4/PT4» » 

COMPUTE PINAL SIGMA*»ZERO ESTIMATES 

SIGNAC I t*VKOB«G4tNSQC I IMO«*«LOGI OCPCI ) >-IO**ALOGI 0< AC I > I ♦ 
•PR9®TC tl«Z«C II 

40 CONTINUE 

tPC IBUG.FO.O) GO T9 A? 

DUNP OS BUG VALUPS TO PRINT SET NBR Z 

WRI' r EC9«41 »TINP v PSUBT*PC4«PTA v PSPeC*R*A*BPAM t AN6TtPR0PT* 
•POO® • DPOOP 

41 PORMATCIH •?CI2**t')*I3«4X» 'CAL PNR*»tF8*2. 

• 9X»«PCAL*4*« .F8«2«/« 

•SX.'PILT PWRt* •BCPS.2I •/« 
l9X»*R(n:S4X«S(P9»2li/» 

•SX**AC II :*t4X*#CPR*ZI*/t 

• 5X* • BCAMC 1 1 S • » 1 X #8CF8* 2 1 •/• 

•SX*«ANGTLC I> :• •BCP4»PI •/• 

•EX ••PR/PTC I » !• •SCP8»2>«/« 

•SXt*FOCI )*• t'*Xt4fF8*2l./. 

•5X**BNDCI >:• «2X«*CFB»?| I 

WRITEI9.I4!) GAf NSO.ZWtSIGMA 
141 FORMAT C 1H *4X • »GRGTC M S • ♦ 1 X • 8CP8« 2) 

•5X» •ROFFC I ) : • *1 X*9CPB*2t«/« 

•5X.*StGMACIi:*«4IF*«2>t 

42 CONTINUE 

BUTLO POINTER SET POR PILLING OUTPUT BUPPSR 
K*M00CMXRCTR*MLI4I 
00 45 J*l.3 
MATRIX! J*Kl* T INPC Jl 
45 CONTINUE 

PUT A/C OATA VALUPS INTO OUTPUT BUFFER 
MA V RI XC4»K)*IPfXC M* 19*1 
MATRIXCS.KI*IPIXCVFL*I0«» 

MATRIXC4«K)«A0ASC3» 

MATRIX! T«K|«A0ASC4| 

MATRIX C B« K I * AOASCS I 

M4 T RIXC9.Kl»IP!XCl«0**ALOGtOCPSU8T| I 
MATQ{xnO*K)>!FIX< 10n*«AL0G14CFNOIZ> » 

POSITION THE 4 SIGMA VALUES IN Th® OUTPUT BUFFER 
on SO 1*1 .8 

I BPTRI I !»IFIX|T ANC ANC.TC | |/ST,3»/THET5PXCELLSP0«5I 
IBK*MXRC t B*!BPT»CI ) 

!F!!BK*LT«0> IBK*0 
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I8K«M00!1BK*MLH>I 

MATRIX! tOM • IBK l«tFt Xl SIGMA! t |9|0« I 
SO CONTINUE 
nout *nout ♦ 1 
S2 mout»mxrctr«lprmt 
C TAKE CARE OF MXRCTR COUNTER OVERFLOW 
fF(N0UT*LT«0» NOUT*MOUT«NL 
MOU r »NOO«MOUr*NLlA! 

IFlNOUT.GT.l I GO TO »0 

DISPLAY FIRST RECORD INFORMATION 

WAITE! 6. 61 »THETA*PDOP*OPDOP*R*A*IBPTR 

61 FORMAT! 1 HO* • RCO FI PARAMETERS* • /• t OX •' ANGLES (OEGt S •*B!SX*F6.1» 
!/• 10X**DORPLER FREOS: • *6!SX*F6»t > • 

*/♦ I 9X* * BAND WIDTHS: • »9(5X,E6.1 I*/ 

* t ok ••ranges: • •ak.mex.fe.i »./, 

6I0X* * AREAS : • *AX*8!SX*F6«t »•/, 

6 1 OX ♦• POINTER VALS: ***61 SX • IF )/ > 

WRITE!6*62» ITHETA 

62 FORMAT(SX**TIMf *,6 X**acT *,IX«*VEL *. 2X *• PITCH* . 3X* • ROLL **3X 

• *DR|FT* * 7X • *PT* • 3X • *RKOI 2* • 3X • *SI GMAI II ♦ I *1 *8*./* 

S62X*S! 12* *OEG* *2X1 *//l 

FORMAT THE OUTPUT FOR PRINTER SET NBP 1 
TO OO 6S I *3* 1 0 

TRIX 1 1*2 1 pFLOATI MATRIX C l *MOUT | )/! 0* 

66 CONTINUE 

OUTPUT THE PRINTEO OATA SET NO© I 

WRI TE« 6 #00M MATRIX! K,MOUT)*X*l *2) .ITR I X! II • I *1 • 1 6 1 

60 FORMAT ! t X* 2! 12** t • I *F4# 1 *1 X* 16!F6« I • I XI I 

OUTPUT THE PUNCHED OATA SET *6R | 

WRITE! T # 81 I(MATRIX!K*M0UT»,K*| *161 

61 FORMAT! 2t 2 • I 3*219*314*1 0191 

8ACX-FILL THE BUFFER LINES THAT HAVE BEEN OUTPUT 

00 oo 3*1*16 

MATRIX! J*MOU T I **999 
90 CONTINUE 

mxpctr*mxrctrm 
C TAKE CARF OF REGISTER OVFR»FLCW 
tP!MXRCTR*EQ«0 I MXRCTR«K#I 
t F| NOUT«GE« I END I GO *0 T 
GO TO 29 
777 CONTINUE 

C PURGE THE BUFFER BEFORE STOPPING 
ICLR»ICLRT1 

IF! ICLR#GT*LPRNT| GO to TTR 
K*MOO ! MXRCTR *ML I ♦ I 


c 

c 

c 


60 TO S2 

77# ««tTet6*'’00l NOU^T!** 

700 FORMAT! t HO tSX»*NRR OF RfCOROS PRO CCS SCO ■ St?,5K. 
• •LAST ACCORD Tim? a • .?< t2 . • : • » ,T 31 
STOP 

#00 WRtTFIB.SOII MLN.ML 

#0t FORMAT! 1H .****> MATRIX BUFFER OVERFLOW <****«16. 
•• LINES NECOEO* • • 1 #• * LINES AVAILABLE* •/ > 

60 TO 77T 
END 

SUBROUTINE PWB|IL.IR*AFT,PSD» 

THIS ROUTIN' USES LEFT 0 P 16HT POINTERS TO CALCULATE 
THE MFAN VALUE OF POWER IN THE FREO CELL specified 


DIMENSION AFT! 10241 
TEMP»0» 

DO 10 I *IL*t R 
TENP'TENP'AF^I | |4*2 
ID CONTINUE 
PSOaTfMP 
RETURN 




END 

SUBROUTINE I NPUT ( PO I N T S • T I 
INTE6ER POINTS 

INTEGER*? t? me C SI • ADA SI • | » SINE I 20481 *COSlNF( 20481 
COMMON/I T APE /▼ l ME • ADAS . SIN*. COSINE 
t9T*aPOlNTR/IO?4 
60 TO II0.201.IPTS 

ID RFA0C8.IS.FNn«T^ ) T\ ** , AOAS . I SI NEC J! • J«I .POINTS I • 

41 COSINE! K | .«■! .POINTS! 

IS FORMA T ( SA2. 1 6! 84 A2 1 • I Af 44 A2 1 I 
GO TO 30 

20 READ! 8 • 2D **ND*?0 I t|mp, *OAS« I SINE! J1 » J« I .POINTS 
*! COSINE!* I .K«l , ©HINTS ! 

2S FORMAT! 9 42.321 84 A? I .321 64 A2 I I 
30 HTbTIMFU » 

MTaTIMEI 21 
STsTIMFI 3» 

T*3800»*MT*A0. *M*>ST/|0. 


RETURN 

TD MRITE<6,T5> 

7S FORM A T 1 1 H . 1 OX . * ENO OF F IL e ON T*Pf gNlT ****•./) 


STOP 


END 

3UDR0UTINE SCALXlPOlNTS.SINE .COSINE.FRS! 

integer POINTS 

INTEGER*2 SIN*! POINTS I .COSINE! POINTS » 

complex frsipointsi .cmplx 

COMMON /RAT/FINCR .FCAL.FSANP. INC# 

COMMON/F ILT/I TyPE .XP0L2 


on on on on 
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C BYPASS AVERAGING WHEN 8*4GHZ DATA 
GO TO (4 *4 *4 .6 1 * I TYPE 
A CONTINUE 

C AVERAGE OUT THE OC COMPONENT 
SUMC*0. 

SUMS*0* 

00 5 M*1 .POINTS 
ISIN*SINF(M) 

IC0$*C0SINF< Ml 
SUMS*SUMS*FLOAT< ISINI 
SUMC«SUNCTELOAT< ICOSI 
8 CONTINUE 

I SAVG«IEtX(SUMS/ELOAT( POINTS)! 

ICAVG*IFIX( SUMC/FLOAT! POINTS) ) 

8 CONTINUE 

C BUILD <XtJV) BUFFER FOR OFT 
00 10 K*l .POINTS 
1SIN*S1NE!K) 
tCOSsCOStNECK) 

C ZERO CH2 WHEN 0*4GHZ DATA 
GO TO < 8. 8.8.7 ) • ITYPP 
T IC0S«0 
t SAVG*0 
!CAVG*0 
B CONTINUE 

IS1N*I SIN-ISAVG 
IC9S«1C0S-ICAVG 

PRS(K)*CMPLX<FLOAT(ISIN/16)/2048*.FLOAT| TC0S/l6)/2048* > 

10 continue 

RETURN 
EN0 

FUNCTION GAMMA(FREO) 

WATER ROLLOFF FILTER FOR CALCULATING ZW!I) 

FREO EXPECTED IN HZ 
COMMON /FILT/I TYPE. KPOLZ 
P*PREQ/1000* 

GO TO (10*20. 30*40). ITYPE 
THIS FUNCTION IT FOR 13.3GHZ 
>> OATA TAKEN AFTER IJAN78 «< 

1 0 GAMMA* 1 5. 8532-F* < T. TABS l-F* < I . 5464-FM . 1 0Q227-* 000945247*F ) > > ♦ 
4(4*621 35-* 1721 3R/FI/F 

15 GAMMA * -GAMMA 

16 RETURN 


THIS FUNCTION IS FOR 1*6 GHZ 

20 GAMMA*7.7|0938-FM 2 1*64141-1 A. RA1A1*F)4< 1* 1 1 1 81 6-0* 01 86787/F) /F 
GO TO 15 

THIS FUNCTION IS FOR 4*75GHZ 

30 GAMMA* l 8* 5575— F4( 1 6* 5Q8 t — F* ( 8*640R7«f* ( 3*44443— 0*65321 3*F ) ) )- 


i 
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K4t2SS4Q/Ft 

go ro is 

▼HIS FUNCTION IS FOR 0*4GHZ 
40 GAMMA* 0*0 
GO TO 16 
END 

SUBROUTINE PLOTIX.V.N) 

INTEGER 0 
INTEGER QA.OB.OC 
INTEGER X.V 

30 FORMAT! • •,IS,SX,IS,SX«1X,103A|) 

30 FORMAT I • • ,2E1 3, S, 1 X, 1 03AI » 

O I MENS t ON Of 1001 «X(N)»Y(N) 

DATA QA,08.0C/2HII,2H**,2H— / 

00 1 1*1.100 
Q( 1 1 =OC 
1 CONTINUE 
VMAX*V«1 ) 

VMINsVll ) 

DO 3 1*2, N 

IF( YMAX.LT, Y( II) YMAX*Y( I ) 

IF(YMIN,GT.Y( I I) YMIN*V(1) 

3 CONTINUE 

SCALE=YMAX-VMTN 
I F< SCALE* EQ. 0 I GO TO 16 
OO 2 1*1 .N 

NUM*< < V( I )-YMTNI*RA,0)/SCALEM 
2 WRITE (6*30) XI I ) . Y( t ) »CA , ( 0( 1 1 ) • 1 1*1 .NUM) ,OB 
RETURN 

16 WRITE 16,33) YMAX 

33 FORMAT (• DEPENDENT VARIABLE IS CONSTANT *,D13,5 ) 

RETURN 
END 

SUBROUTINE OFT | X , M , t NCR ) 

C A FAST FOURIER TRANSFORM THAT IMPLEMENTS THE 

C "DECIMATION-IN-FREOUENCV" algorithm 

C 

C CALLING PARAMTEOS: 

C x * complex array of dimension n that CONTAINS INITIALLY 

C the input SEQUENCE XCN) AND FINALLY CONTAINS THE 

c TRANSFORM X ( K ) • THE OUANTITY M is AN INTEGER, M*L0G2(N), 

C WHERF N IS the NUMBER OF TIME VARYING DATA POINTS 

C l NCR * SIZE OF THE ARRAY MOLDING THE time FUNCTION 

C 

COMPLEX X( INCR),U.W.T,CMPLX 

N*2«*M 

NV2*N/2 

NM!*N»1 

J*1 



00 7 1*1 *NM| 

IF(t.GE.J) GO TO 5 
Tjs*.| J) 

*1 J)*X<t ) 

X<I)*T 
S K«NV2 

ft tF(K«GE* J ) GO TO 7 
J*J*K 
K*K/2 
GO TO 6 
7 J*J4K 

■5 PI*3«14S02ft5 
00 20 L*1.M 
LS*t**L 
LEl*LE/2 
U*CMPLX( 1 »0»0«0 I 

M*CMPLX<COS< PI/FLOAT <LP| >> .•StNCP!/FLOATU. e l » » ) 

00 20 Jsl.LF! 

OO 10 I*J.N,LF 
1P*ItLE 1 
TaXf IP)7U 
X(IPI*X(t)«T 
10 Xftl*XCII*T 
20 U*U*V 
RETURN 
END 

SUBROUTINE CARFAlAREA.BEAM.MFlGT.ANGLO.VELOe. ROLL. DRIFT. ICON) 

REAL KONE.KTWO 

COMMON /KAREA/KONE.KTmO 

TERM|*2« *T AN< BFAM/2. I*HE IGT**2/C0S< ROLL > /CClSI ANGLO > 

TEMPI =KTwO*VFLOC**2 

ATO*TANlORlFT) 

TNR*T AN( ROLL I 
COR*COi»< DRIFT) 

TEMPI a|-AT0nTNR*SQRTCTNR**2*(TCMPi«i # 0)+ 

• TEMPI *CDR**2-1 • 0 ) ) /< t • 0-TEMPI *CDR»*2 ) 

TERM2aTEMPl 

TEMPI =KONE*V»LOC**2 

TEMPI s«-AT0*TNP*SQRT<TNR**2* I TEMP 1»1*0 »♦ 

• TEMP l* COR** 2*1 • 0 I ) /< 1*0* TEMPI *CDR**2) 

IF| ICON. LT. 01 TEMPl^O. 

^RFA*TERM1*(TERM2-TEMPI I 

ARSA«ABS< AREA ) 

RETURN 

END 
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C GAIN SOB- PR OGB AH FOR 13.3GHZ VV DATA 
BBAL FUNCTION 6AZN (ANTING, PITCH) 

REAL ANTNNA, ANGLE 
REAL ANT ANG, PITCH* 

DIRE IS ION ANTNNA (71) 

GAIN TALOBS FOR 13.3 RHC ANTENNA# VERTICAL POLARIZATION 
»» TABLE OP-DATED 16JAR1979 ...«< 

DATA ARTNNA(OOI) #AHTNHA(002) #AHTNNA(003) #ANTNRA (004) #ANTNNA (005) # 

* ANTRNA (006) #ARTNNA (007) # ANTNNA (008) #AHTRNA (009) #ANTNRA (010) 
4/17.00,17.00,18.30,19.50,21.00, 

* 22. 50 #23. 80 #24. 70,26.00 #27. 50/ 

DATA ANTRNA (01 1)#ANTNNA (012)# ANTNNA (013) , ANTNNA (014) , ANTNNA (015) # 

* ANTNNA (016) , ANTNNA (017) , ANTNNA (018) # ANTRNA (019) , ANTNNA (020) 
4/28.80, 30.00,30.50, 30. 70#31.00, 

* 31.00, 31.50, 32.00,33. 00,34.00/ 

DATA ANTRNA (021) , ANTNNA (022) , ANTNNA(023) , ANTNNA (024) , ANTRNA (025) , 

* ANTRNA (026) # ANTNNA (027) # ANTRNA (028) , ANTRNA (029) , ANTRNA (030) 
4/34.40,34. 80# 35. 30, 35. 50,35.50# 

* 35.50,35.30,35.00,34.00,33.20/ 

DATA ARTBNA(031) # ANTNNA (032) , ANTRNA (033) .ANTNNA (034) , ANTNNA (035) , 

* ANTNNA (036) .ANTNNA (037) .ANTRNA (038) .ANTNNA (039) .ANTRNA (040) 
4/33.00,33.50,34.00,34. 50,35.00, 

4 35.50,35.00,34.50,34.00,33.50/ 

DATA ANTNNA (041) .ANTRNA (042) , ANTNNA (043) .ANTNNA (044) .ANTNNA (045) , 
4 ANTRNA (046) .ANTNNA (047) .ANTNNA (048) .ANTNNA (049) .ANTRNA (050) 

4/33.00,33.20,34.00,35.00,35.30, 

4 35. 50, 35. 50, 35. 50, 35. 30, 34. .90/ 

DATA ANTNNA (051) .ANTNNA (052) .ANTNNA (053) .ANTNNA (054) , ANTNNA (055) , 

* AHTNNA(056) .ANTRNA (057) .ANTNNA (058) .ANTRNA (059) .ANTNNA (060) 
*/34. 40, 34. 00, 33. 00, 32. 00,31.50, 

4 31.00,31.00,30.70,30.50,30.00/ 

DATA ANTNNA (061) .ANTRNA (062) , ANTNNA (063) .ANTNNA (064), ANTRNA (065) , 
« ANTNNA (066) .ANTNNA (067) .ANTRNA (068) .ANTNNA (069) .ANTNNA (070) 
4/28.80,27.50,26.00,24.70,23.80, 

* 22.50,21.00,19.50,18.30,17.00/ 

DATA ANTNRA(071) /17.00/ 

ANGLE IS THE DIFFERENCE OF THE ANTENNA ANGLE AND PITCH 
ANTING IS THE ANTENNA ANGLE BITH RESPECT TO THE AIRCRAFT 
PITCH IS THE AIRCRAFT PITCH ANGLE 

ALL ANGLES ARE EXPECTED TO BE EXPRESSED IN RADIANS 

THE RANGE OF VALUES OF ANGLE IS: -70.9 TO 70.9 DEGREES. 

IF AN ANGLE FALLS OUT OF THIS RANGE, THE CLOSEST AHGLE 
VALUE HILL BE RETURNED 

ANGLE 3 ANT ANG-PITCH 
INDEX- IFIX ( (ARGLB*57. 3470. )/2.) 41 
IF (INDEX. LE.O) GO TO 30 


n no 


XPCZIon.GS.71) GO TO 10. 

ARGL2* FLOAT (IHDBX-1) *2.-70. 

GAXR«AHf RHA (IHOBX) ♦ (AHTHHA (XRDBX* 1) -AHTHHA (XIDBX) ) • 
S | (ARGLB*57. 3- A HOLS) /2 • 0) 

BBTOBH 

10 GAXR«ARTBRA (71) 

BBTOBH 

30 GAZR* AHTHH A ( 1 ) 

BBTOBH 

BHD 


BBAL POHCTXOH GAXH (AHTARG, PITCH) 

BBAL AHTHHA, AHGLB, AHTARG, PITCH 
DXRBHSIOH AHTHH A (71) 

GAIH fALOBS FOB 1.6 HOC A HTENH A, VERTICAL POLARIZATXOH 

DATA AHTHH A (001) ,AHTHHA (002) , AHTRIA (003) , AHTHHA (004) , AHTHHA (005) , 

* AHTHHA (006) , AHTHHA (007) , AHTHHA (008) , AHTHHA (009) , AHTHHA (010) / 

* 18.70, 15.40, 15.70, 16.40, 17.40, 

* 17.90, 18.70, 19.40, 19.90, 20.60/ 

DATA AHTHHA (Oil) , AHTHHA (012) , AHTHHA (013) ,ANTHHA (014) , AHTHHA (015) 

* AHTHHA (016) , AHTHHA (017) , AHTHHA (018) , AHTHHA (019) , AHTHHA (020) 

* 21.60, 21.90, 21.90, 22.40, 22.70, 

* 22.90, 22.60, 22.90, 22.40, 22.20/ 

DATA AHTHHA (021) ,AHTHHA(022) , AHTHHA (023) , AHTHHA (024) , AHTHHA (025) 

* AHTHHA (026) , AHTHHA (027) , AHTHHA (028) , AHTHHA (029) , AHTHHA (030) 

* 21.90, 21.40, 21.40, 21.10, 20.90, 

* 20.70, 20.40, 20.40, 20.10, 20.70/ 

DATA AHTHHA(031) , AHTHHA (032) , AHTHHA (033) , AHTHHA (034) , AHTHHA (035) 

* AHTCHA (036) , AHTHHA (037) , AHTHHA (038) , AHTHHA (039) , AHTHHA (040) 

* 20.70, 21.40, 21.40, 21.70, 22.10, 

* 22.20, 22.20, 22.10, 21.90, 21.40/ 

DATA AHTHRA(O'I) , AHTHHA (042) , ARTHHA(043) , AHTHHA (044), AHTHHA (045) 

* AHTHHA (046) , AHTHHA (047) , AHTHHA (048) , AHTHNA (049) , AHTHHA (050) 

* 21.20, 20.90, 20.70, 20.10, 19.70, 

* 19.20, 18.70, 18.10, 17.60, 17.60/ 

DATA AHTHHA(OSI) , AHTHHA (052) ,ARTHHA(053) , AHTHHA (054) , AHTHHA (055) , 

* AHTHHA (056) , AHTHH A (057) , AHTHHA (058) , ANTHHA (059) , AHTHHA (060) / 

* 17.40, 17.10, 17.20, 17.10, 16.90, 

* 16.90, 16.90, 16.70, 16.40, 15.90/ 

DATA AHTHHA (061) , AHTHHA (062) , AHTHHA (063) , AHTHHA (064) , AHTHHA (065) , 

* AHTHHA (066) , AHTHHA (067) , AHTHHA (068) , AHTHHA (069) , AHTHHA (070) / 


15.40, 14.90, 14.70, 

14.20, 

13.90, 

13.20, 12.70, 11.40, 

10. 40, 

9.40/ 

DATA AHTHHA (71) /8. 90/ 
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THIS SUBROUTINE RETURNS A GAIN VALUE FOB IB ARGLB 

AH6LB IS TBB DIPPBBBNCB OP TRB ARTBBNA ARGLB ABO THE PITCH 
AHTAHG IS TRB ABTBBBA ARGLB BITR RBSPBCT TO TRB AIRCRAFT 
PITCH IS TRB AIRCRAFT PITCH ARGLB 

ALL ARGLBS ARB BXPBCTBD TO BB EXPRESSED IB RADIABS . 

THE RABGB OP VALUES FOR ARGLB IS: -70.9 TO 70.9 DBSREES 

IF AR ARGLB FALLS OUT OF THIS BANGS, TRB CLOSEST GAIN 
VALUE HILL BB RBTUBRBD 

ARGLB * ARTARG-PITCB 
INDEX* IPIX ( (ARGLF*57. 3 ♦70. ) /2. ) ♦ 1 
IF (INDEX. LB. 0) GO TO 30 
IF(IHDBX.GB.71| GO TO 10 
ANGLZ*FLOAT(INDBX-1) *2.-70. 

GAIR-AHTRHA (INDBX) ♦ (ANTNNA (INDEX* 1) -ANTNRA (INDEX) ) *ABS 
S ( (ANGLE*57. 3-ARGLZ) /2.0) 

BVTtlUM 

10 GAIR*ANTNNA (71) 

RETURN 

30 GAIH*ANTNNA(1) 

RETURN 

END 


REAL FUNCTION GAIN (ANTING* PITCH) 
REAL ANTNRA* ARGLB* ANTING .PITCH 
DINENSION ANTNNA (71) 


GAIN VALUES FOR 1.60 KHC ANTENNA* CROSS POLARIZATION (V-H) 


I 

i 

$ 

i 

i 

i 


DATA ANTNNA (001) * ANTNNA (002) , ANTNRA (003) .ANTNNA (004) .ANTNRA (005). 

* ANTNNA (006) .ANTNNA (007) * ANTNRA (008) .ANTNNA (009) .ANTNNA (010) / 

* 21.40* 22.20* 22.70. 23.20. 23.90. 

* 24.10. 24.60* 24.90. 25.10. 25.90/ 

DATA ANTNNA (0 11) .ANTNNA (012) *ANTNNA(013) * ANTNNA (014) .ANTNNA (015) . 

* ANTNNA (016) .ANTNNA (017) .ANTNNA (018) .ANTNNA (019) * ANTNNA (020)/ 

* 26.10. 26.20. 26.20. 26.40. 26.20* 

* 26.20. 25.60. 25.60. 24.60, 24.20/ 

DATA ANTNNA (021) .ANTNRA (022) .ANTNNA (023) .ANTNNA (024) .ANTNNA (025) 

* ANTNNA (026) .ANTNNA (027) .ANTNNA (028) .ANTNNA (029) .ANTNNA (030) 

* 23.60. 22.90. 22.20. 21.60, 20.90* 

* 20.10. 20.10, 19.70. 19.40. 19.70/ 

DATA ANTNNA (031) .ANTNNA (032) .ANTNNA (033) .ANTNNA (034) , ANTNNA (035) . 

* ANTNNA (036) .ANTNNA (037) ,ANTNNA(038) .ANTNNA (039) , ANTNNA (040) / 

* 19.90* 20.10* 20.90, 20.90* 21.40, 

* 21.70, 22.20, 22.20* 22.20, 21.90/ 

DATA ANTNNA (04 1) .ANTNRA (042) ,ARTNNA(043) .ANTNNA (044) * ANTNRA (045) , 



x 
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* AITIIA (046), AITIIA (047) .AITIIA (048) .AITIIA (049) .AITIIA (050) / 

* 21.90, 20.90, 20.70, 20.10, 19.70, 

* 19.70, 20.20, 20.40, 20.70, 20.90/ 

Dm A If ill (051) .AITIIA (052) ,A0T10A(053) ,A«TIIA (054) , AITIIA (055) 

* AITIIA (056) , AITIIA (057) , AITIIA(058) , AITIIA (059) , AITIIA (060) 

* 21.20, 21.40, 21.40, 21.40, 21.20, 

* V '60, 20.90, 20.40, 19.90, 19.10/ 

DATA AITIIA (061) , AITIIA (062) , AITIIA (063) .AITIIA (064) .AITIIA (065) 

* AITIIA (066) .AITIIA (067), AITIIA (068) .AITIIA (069) .AITIIA (070) 

* 18.40, 18.10, 17.70, 17.40, 17.90, 

* 17.60, 17.90, 17.40, 16.60, 16.10/ 

DATA AITIIA (71) /15. 60/ 

THIS SUBBOOTIIB 8BT08IS A GAIB 7AL0B FOB Al AI6LB 

AI8LB IS TIB DIPPBRBICB OP THB AITBIIA AIGLB AID TIB PITCH 
AITAIG IS THB AITBIIA AIGLB IITH BBSPBCT TO TBE AIBC8APT 
PITCH IS THB AIBCBAPT PITCH AIGLB 
ALL AIGLBS ABB BXPBCTBD TO BB BZPBBSSBD II RADIUS 

THB BAIGB OP VALOBS POR AIGLB IS: *70.9 TO 70.9 DBGBBBS 

IP Al AIGLB PALLS OOT OP THIS BAIGB, THB CLOSEST GAII 
VALUE HILL BB RBTOBIBD 

AIGLB * AITAIG-PITCH 
IIDBI* IFIX ( ( AIGLB*57. 3 ♦TO. ) /2. ) ♦ 1 
IP (IIDBI . LB. 0) GO TO 30 
IP (IIDBI. GB.71) GO TO 10 
AIGLZ* FLOAT (IIDEX-1) *2.-70. 

G All* AITIIA (IIDBX) ♦ (AITIIA (IIDBX+ 1) -AITIIA (IIDBI) ) «ABS 
S ( (ABGLB*57. 3-ABGLB) /2 . 0) 

RETORI 

10 GAII* AITIIA (71) 

RBTORI 

30 GAII*AITBHA(1) 

RBTORI 
BID 


** AITBIIA GAII FOICTIOI POR 1.6GHZ HH DATA 
BBAL FOICTIOI GAII (AITAIG, PITCH) 

RBAL AITIIA, AIGLB, AITAIG, PITCH 
DIHBHSIOH AITIIA (71) 

GAII PALOBS POR 1.6 KHC AITBIIA, HORIZOITAL POLABIZATIOI 

DATA AITIIA (001), AITIIA (002) , AITIIA (003) .AITIIA (004) , AITIIA (005) , 

* AITIIA (006) .AITIIA (007) ,AITIRA(008) .AITIIA (009) .AITIIA (010) / 

* 23.10, 23.70, 24.20, 24.40, 24.40, 


«(\ 
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• 24.40, 24.60, 24.60, 24.60, 24.90/ 

DATA AVTMA(OII) ,ANTIVA(012) , AVTVVA(013) ,AITVIA (014) ,AKTR«A (015) , 
4 AVTIIA (0 1 6) , AITIf A (0 1 7) , AVTVR1 (018), AITOVA (0 1 9) , A IT! VA (020) / 

• 24.60, 24.40, 24.40, 24.40, 23.90, 

• 23.70, 23.10, 22.90, 22.40, 21.90/ 

DATA AVTIIA (021) , AITVIA (022) , AVTIIA (023) , AVTIIA (024) ,AITIIA (025) , 
4 AITVIA (026) , AVTIIA (027), AVTIIA (028) , AVTIIA (029) ,AITIIA (030)/ 
4 21.40, 20.90, 20.40, 19.90, 19.10, 

4 18.90, 18.60, 18.40, 18.40, 18.40/ 

DATA AVTIIA (031) , AITVIA (032) , AVTIIA (033) , AITVIA (034) , AVTIIA (035) , 
4 AITVIA (036) ,AVTVIA (037), AVTIIA (038) , AITVIA (039) , AITVIA (040) / 

• 18.90, 18.90, 19.70, 19.70, 19.70, 

4 19.70, 19.70, 18.90, 18.70, 18.10/ 

DATA AiriIA(041) .AITVIA (0*2) .AITVIA (043) , AVTIIA (044) , AVTIIA (045) 

• AITVIA (046) , AVTIIA (047), AVTIIA (048), AITVIA (049), AVTIIA (050) 

• 17.60, 17.70, 17.90, 18.40, 18.20, 

• 17.90, 17.90, 18.10, 19.40, 19.90/ 

DATA AVTIIA (051) , AITVIA (052) , AITVIA (053) .AITVIA (054) .AITVIA (055) 

4 AITVIA (056) .AITVIA (057) .AVTIIA (058) .AITVIA (059) .AITVIA (060) 

• 20.70. 21.20. 20.90. 21.20. 21.90. 

• 23.10. 23.40. 23.40. 22.90. 22.20/ 

DATA AV?VIA(061) .AVTIIA (062) .AITIIA(063) .AITVIA (064) .AITVIA (065) 

4 AITVIA (066) .AITVIA (067) ,AVTIIA(068) .AITVIA (069) .AITVIA (070) 

4 21.90, 22.10, 22.20, 22.20, 22.40, 

4 22.10, 21.90, 21.60, 21.90, 22.10/ 

DATA AITVIA (71) /22. 10/ 


THIS SUBROUTINE RBTORVS A GAIV VALUE FOB AN AISLE 


AV6LE IS THE DIFFERENCE cF TBB AVTBVVA ANGLE AND THE PITCH 
AVTAVG IS THE AITEVVA ANGLE WITH RESPECT TO THE AIRCRAFT 
PITCH IS THE AIRCRAFT PITCH ANGLE 
ALL ANGLES ARB EXPECTED TO BE EXPRESSED IN RADIANS 


I 

f 

I 

i 


» 
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i 
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THE BARGE OF VALUES FOR ANGLE IS: -70.9 TO 70.9 DEGREES 

IF AV ANGLE FALLS OUT OF THIS RANGE, THE CLOSEST GAIV 
VALUE HILL BE RETURNED 

1 

AVGLE « AVTAVG-PITCH 
IRDEX>IFIX ( (ARGLE457. 3*70. ) /2. ) M 
IF (INDEX. LE.O) GO TO 30 
IF(IRDBX.GE.71) GO TO 10 
AVGLZ*FLOAT ( INDEX- 1 ) *2. -70. 

GAIV*ARTVVA (INDEX) ♦ (AHTRRA (I NDEX+ 1) -AVTIIA (INDEX) ) *ABS 
$ ( (ANGLE457. 3-ARGLZ) /2 .0) 

RETURN 

10 GAIH*ARTHHA (71) 

RETURN 

30 GAIV>ARTVVA(1) 

RETURN 
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C M ARTBHHA GAIB FOBCTXOB FOR 1.6GHS HV DATA 
RBAL FOBCTXOB OAXB (ART AUG, PITCH) 

SEAL ABTBBA, ABGLB.AHTAHO, PITCH 
DXHBRSXOB ABTBBA (71) 

6AXB TALOBS FOR 1.60 KflC ABTEHRA, CROSS PLOARXZATXOB (H-7) 

DATA ABTBBA (001), ABTBBA (002) , ABTBBA (003) , ABTBBA (00O) # ABTBBA (005) , 

* ABTBBA (006) , ABTBBA (007) , ABTBBA (008) , ABTBBA (009) , ABTBBA (010)/ 

* 16.00, 17.20, 17.00, 17.70, 18 ,20, 

* 18.00, 18.70, 19.00, 19.60, 19.90/ 

DATA ABTBBA (0 1 1) , ABTBBA (012) , ABTBBA (0 1 3) ,AHTBtTA (010), ABTBBA (015) , 

* ABTBBA (016) , ABTBBA (01 7) , ABTBBA (016) , ABTBBA (019) , ABTBBA (020) / 

* 20.10, 20.20, 20.00, 20.00, 20.00, 

* 20.00, 20.60, 20.00, 20.20, 20.10/ 

DATA ABTBBA (021) , ABTBBA (022) , ABTBBA (023) , ABTBBA (020) , ABTBBA (025) , 

* ABTBBA (026), ABTBBA (027) , ABTBBA (028) , ABTBBA (029) , ABTBBA (030) / 

* 20.20, 19.90, 19.90, 19.90, 19.00, 

* 19.00, 19.20, 19.00, 19.00, 19.90/ 

DATA ABTBBA (031) , ABTBBA (032) , ABTBBA (033) , ABTBBA (030) , ABTBBA (035) , 

* ABTBBA (036) , ABTBBA (037) , ABTBBA (038) , ABTBBA (039) , ABTBBA (000)/ 

* 20.20, 20.00, 20.70, 20.70, 20.00, 

* 20.20, 19.70, 18.90, 18.00, 17.60/ 

DATA ABTBBA (001), ABTBBA (002) , ABTBBA (003) , ABTBBA (000) , ABTBBA (005) , 

* ABTBBA (006) , ABTBBA (007), ABTBBA (008) , ABTBBA (009) , ABTBBA (050)/ 

* 17.00, 17.90, 17.90, 18.00, 18.20, 

* 17.00, 16.00, 16.40, 16.90, 16.90/ 

DATA ABTBBA (051) , ABTBBA (052) , ABTBBA (053) , ABTBBA (050) , ABTBBA (055) , 

* ABTBBA (056) , ABTBBA (057) , ABTBBA (058) .{‘BTRBA (059) .ABTBBA (060) / 

* 16.90, 16.90, 16.70, 17.20, 18.20, 

« 19.00, 19.00, 19.70, 19.00, 18.90/ 

DATA ABTBBA (061) , ABTBBA (062) , ABTBBA (063) , ABTBBA (060) , ABTBBA (065) , 

* ABTBBA (066), ABTBBA (067), ABTBBA(068) .ABTBBA (069) .ABTBBA (070) / 

* 18.90, 19.10, 19.20, 18.90, 18.00, 

* 17.00, 16.70, 16.10, 15.90, 15.60/ 

DATA ABTBBA (71) /15. 00/ 

THIS SOBROOTIBB RBTORBS A GAIB VALOS FOR AR AB6LE 

ABGLB IS TBB DIFFBRBBCB OP THE ABTEBBA ABGLB ABD TBB PITCH 
AHTABG IS TBB ARTBHHA ABGLB VITH RESPECT TO THE AIRCRAFT 
PITCH IS THB AIRCRAFT PITCH ABGLB 
ALL ABGLES ARE EXPBCTBD TO BE EXPRESSED IB RADIABS 

THB RARGB OF VALDES FOR ABGLB IS: *70.9 TO 70.9 DEGRESS 

IF AH ABGLB FALLS OOT OF THIS RAHGE, THB CLOSEST GAIB 



on n on 


VALUE BILL BE RETURNED 


ANGLE • ANTANG-PITCH 
XNDSX*IPXX ( (ANGLE+57. 3»70.)/2.) *1 
IP (INDEX. LE.O) GO TO 30 
IP (INDEX. OB. 71) GO TO 10 
ANGLE* PLOAT (INDEX* 1 ) *2. -70. 

GAXN-ANTNNA (INDEX) ♦ (ANTRNA (INDEX* 1) -ANTNNA (INDEX) ) *A8S 
S ((ARGLB*57.3-ANGLZ)/2.0) 

RETURN 

10 GAIR*ANTNNA(71) 

RETURN 

30 GAIN*ANTRNA(1) 

RETURN 

END 


GAIN VALUES POR .400 NEC ANTENNA, VERTICAL POLARIZATION 
REAL FUNCTION GAIN (ANTANG# PITCH) 

REAL ANTNNA# ANGLE# ANTANG# PITCH 
DIHENSION ANTRNA (71) 


DATA ANTNNA (001) 

* ANTNNA (006) 

* 10.40# 12. 

* 19.40# 20. 

DATA ANTNNA (Oil) 

* ANTNNA (016) 

* 23.80# 24. 

* 24.40# 23. 

DATA ANTNNA(021) 

* ANTNNA (026) 

* 19.40# 18. 

* 11.40# 9. 

DATA ANTNNA (031) 

* ANTNNA (036) 

* 4.80# 4. 

* 2.80# 1. 

DATA ANTNNA(041) 

• ANTNNA (046) 

* -1.60# -2. 

* -4.60# -4. 

DATA ANTNNA(OSI) 

* ANTNNA (056) 

* -1.60# -1. 

* - 0 . 20 # - 0 . 

DATA ANTNNA(061) 
• ANTNNA (066) 


# ANTNNA (002) 
#ANTNRA (007) 
40# 14.80 

80# 21.80 

# ANTNNA (012) 

#ANTRNA (017) 
40# 24.80 

80# 22.80 
#ANTNNA (022) 
#ANTRHA(027) 
40# 16.80 

80# 7.80 

# ANTNNA (032) 
#ANTNNA (037) 
40# 3.80 

80# 1.40 

#ANTRNA (042) 

# ANTNNA (047) 

20# -3.20 

60# -4.20 

#ANTNNA (052) 
# ANTNNA (057) 
20# -0.60 
60# -0.60 

# ANTNNA(062) 
#ANTNNA (067) 


ANTNNA (003) 
ANTNNA (008) 
16.40# 

, 22.80# 

# ANTNNA (01 3) 
# ANTRNA (018) 

# 24.80# 
22.40# 

ANTNNA (023) 
ANTNNA (028) 
15.40# 

# 6.80# 

# ANTNNA (033) 
# ANTNNA (038) 

# 3.40, 

, 0.40# 

# ANTNNA (043) 
# ANTNNA (048) 

# -4.20# 

# -3.20# 

# ANTNNA (053) 
# ANTNNA (058) 

# -0.60# 

# - 1 . 20 # 

# ANTNNA (063) 
, ANTNNA (068) 


#ANTNNA (004) 
#ANTNNA (009) 
18.40# 
23.40/ 

# ANTNNA (014) 
#ANTNRA (019) 
24.80# 
21.40/ 
#ANTNRA (024) 
#ANTNNA (029) 
13.40# 
5.80/ 

# ANTRNA (034) 
#ANTNNA (039) 
2.80# 
-0.60/ 
#ANTNNA (044) 
# ANTNNA (049) 
-4.60# 
-2.60/ 
#ANTNNA (054) 
# ANTNNA (059) 
—0 . 60 # 
- 1 . 20 / 

# ANTRNA (064) 
# \NTNNA (069) 


# ANTRNA (005) # 
# ANTNNA (010)/ 


, ANTNNA (015) 
#ANTNNA (020) 


# ANTNNA (025) 
#ANTNNA (030) 


# ANTRNA (035) 
, ANTNNA (040) 


# ANTRNA (045) 
#ANTNNA (050) 


# ANTNNA (055) 
# ANTNNA (060) 


# ANTRNA (065) 
, ANTRNA (070) 
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• -2.20, -2.60, -2.60, 

* —0.60, — 0.20, -0.60, 

DATA ARTRRA(71)/-15.60/ 


-0.60, 

- 11 . 20 , 


-5.60, 

-13.60/ 


THIS SOBHOOTIRB RRTORRS A GAIN VALOR FOR AH AHOLB 


AROLB IS THB DXFFRRBRCB OP THB ARTBRRA ARGLB AMD TUB PITCH 
ARTANG IS THB ARTBRRA ARGLB VXTH RBSFBCT TO THB AIRCRAFT 
FITCH IS THB AIRCRAFT PITCH ARGLB 
ALL ARGLBS ARB BXPBCTBD TO BE BXPRBSSBD IR RADIANS 

THB BABGB OF VALOBS FOR ARGLB IS: -70.9 TO 70.9 DBGRBES 

IF AH ARGLB FALLS OOT OF THIS RAHGB, THB CLOSEST GAIH 
VALUE BILL BE RBTORHBD 


ARGLB * ARTAHG-PITCH 
IRDBX* XFXX ( (AHGLB*57. 3*70. ) /2. ) ♦ 1 
IF (IRDBX.LB. 0) GO TO 30 
IF (INDEX. 6B.71) GO TO 10 
ANGLE* FLOAT (IHDBX- 1) *2.-70. 

GAIR-ARTRRA (INDEX) * (ARTRHA (INDEX* 1) -ARTNRA (IRDBX) ) *ABS 
S ( (ARGLE*57. 3-ARGLZ) /2. 0) 

BBTORR 

10 GAIN*ANTNRA(71) 

RBTORR 

30 GAIR*ARTBRA(1) 

RBTORR 

BHD 


GAIH VALOBS FOB .400 KHC ARTBRRA, HORIZONTAL POLARIZATION 
RBAL FUNCTION GAIN (ARTANG, PITCH) 

RBAL ARTRHA, ARGLB, ARTARG, PITCH 
DIHERSXOK ARTNRA (71) 


DATA ARTNRAfOM) .ARTRHA <002) , ARTRHA (003) , ARTRHA (004) , ARTRHA (005) 

* ARTHNA (006) ,ARTHHA (007) .ARTRHA <008) , ARTNRA (009) .ARTRHA (010) 

* 29.40. 29.40. 29.40, 29.40, 29.40, 

* 28.80, 28.80, 27.80, 27.40, 26.80/ 

DATA AHTRRA(OII) , ARTRHA (012) ,AHTRRA (013) .ARTRHA (014) , ARTRHA (015) 

* AHTRRA(016) .ARTNRA (017) , ARTRHA (018) .ARTRHA (019) .ARTRHA (020) 

* 26.40, 25.40, 24.40, 23.40, 21.80, 

* 20.80, 19.40, 17.80, 15.80, 13.80/ 

DATA AHTBR A (021) .ARTNRA (022) , ARTRHA (023) , ARTRHA (024) , ARTRHA (025) , 

* ARTRHA (026), ARTRHA (027) .ARTRHA (028) .ARTRHA (029) .ARTRHA (030) / 

* 11.80, 9.80, 7.80, 5.80, 3.80, 

* 2.40, 0.40, -0.60, -1.60, -2.20/ 

DATA ANTNRA(031) .ARTRHA (032) .ARTRHA (033) .ARTRHA (034) , ARTRHA (035) , 


ob 




* ARTRRA <036) <037) , ARTRRA (038) , ARTRRA <039) , ARTRRA (090) / 

* -2.00, -2.00, -2.60, -3.20, -3.20, 

* -3.00, -9.20, -9.20, -9.60, -5.60/ 

DATA AITI8A<091) ,ARTRRA (092) , ARTRRA (093) , ARTRRA (099) ,ARf RRA (095) , 

* ARYRH A (096) , ARTRRA (097) , ARTRRA (098) , ART RRA (099) , ARTRRA (050) / 

* -6.20, -7.20, -8.20, -8.60, -9.20, 

* -9.60, -9.60, -9.60, -9.20, -8.60/ 

OATA ARTRRA (051) , ARTRRA (052) , ARTRRA (053) , ARTRRA (059) , ARTRRA (055) 

* ARTRRA (056), ARTRRA (057) ,ARTRHA(058) , ARTRRA (059) , ARTRRA (060) 

* -8.20, -7.60, -6.60, -6.20, -5.60, 

* —5.20, —9.60, —9.60, —9.20, —3.60/ 

OATA ARTRRA (061), ARTRRA (062) , ARTRRA <0li3) ,ARTRRA (069) , ARTRRA (065) 

* ARTRRA (066) , ARTRRA (067) , ARTRRA (068) , ARTRRA (069) , ARTRRA (070) 

* -3.60, -3.60, -3.60, -3.60, -9.20, 

* -9.20, -9.60, -5.20, -5.20, -6.20/ 

OATA ARTRRA (71) /-6.60/ 

THIS SQBBOOTXRB RBTORRS A 6AIR 7AL0B FOR AR AR6LB 

AHGLB IS THB DIFFER BRCB OF TUB ARTBRRA AR6LB ARO TUB FITCH 
ARTARO IS TRB ARTBRRA AROtB RITR RBSFBCT TO THB AIRCRAFT 
FXTCB IS TBB AIRCRAFT FITCH AROLB 
ALL AR6LBS ARB BIFECTEO TO BE BIPRBSSBD IR RADZARS 

TBB RAR6B OF VALOBS FOR AR6LE IS: -70.9 TO 70.9 0B6RBBS 

IF AH AR6LB FALLS OOT OF THIS RAROB, TBB CLOSEST 6AIR 
f ALOE HILL BE RBTORHBD 

AHGLB - ARTARG-PITCH 
IPDBX*IFfX((ARGLB*57» 3 ♦70.) /2.) *1 
IF(IHDBX.LB.O) GO TO 10 
IF (IRDBX.GB.71) GO TO 10 
AHGLS«FLOAT(IHDBX-1) *2.-70. 

GAXH*AHTRHA (IRDBX) ♦ (ARTRRA (IHDEX* 1) -ARTRRA (IRDEX) ) *A8S 
S ( (ARGLB+57. 3-AHGL8) /2.0) 

RBTORR 

10 0AIR-ABTRHA(71) 

RBTORR 

30 GAIH>ARTRRA(1) 

RBTORR 
ERD 


GAIR VALOBS FOR .900 KHC ARTBRRA, CROSS FOLARIZATIOH 
REAL FOHCTIOR GAIR (ARTARG, PITCH) 

RBAL ARTRRA, AHGLB, ARTARG, PITCH 
DIHBRSIOR ARTRRA (71) 
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DATA ARTRRA (001) .ARTRRA (002) .ARTRRA (003) .ARTRRA (000) .ARTRRA (00*) 

• ARTRRA (006) .ARTRRA (007) , ARTRRA (006) , ARTRRA (009) , ANTON A (010) 

• 19.90, 20.90, 22.10, 22.90, 23.90, 

• 29.10, 29.00, 29.80, 2S.10, 25.10/ 

DATA AiTViA(OII) ,A>T8RA (012) , AIVTItlfA (0 1 3) , ARTRRA (019) , ARTRRA (015) 

• ARTIIA(OIO) , AITII9 A (017) ,ARTRHA(01S) ,ARTRRA (019) ,AMTNMA (020) 

• 25.10, 29.90, 29.00, 29.10, 23.30, 

• 22.60, 21.90, 20.30, 10.90, 17.60/ 

DATA ARTRRA(021) ,ARTRRA(022) ,AHTRRA (023) ,ARTRRA (029) , ARTRRA (025) 

• ARTRRA (026) , ARTRRA (027) , ARTRRA (028) .ARTRRA (029) .ARTRRA (030) 

• 15.60. 19.10. 12.30, 10.60, 8.60, 

• 6.90, 5.10, 3.60, 2.60, 1.90/ 

DATA ARTRRA (031) .ARTRRA (032) , ARTRRA (033) .ARTRRA (039), ARTRRA (015) 

• ARTRRA (036) .ARTRRA (037), ARTRRA (038) .ARTRRA (039) .ARTRRA (040) 

• 1.10, 0.90, 0.60, 0.10, -0.20, 

• -0.40, -1.20, -1.40, -2.10, -3.10/ 

OATA ARTRRA (04 1) .ARTRRA (042) , ARTRRA (043) , ARTRRA (044) .ARTRRA (045) 

• ARTRRA (046) .ARTRRA (047) .ARTRRA (048) .ARTRRA (049) .ARTRRA (050) 

• -4.10, -5.10, -5.90, -6.60, -7.10, 

• -7.40, -7.40, -6.90, -6.20, -5.60/ 

DATA ARTRRA (051) .ARTRRA (052) .ARTRRA (053) .ARTRRA (0S4) .ARTRRA (055) 

• ARTRRA (056) .ARTRRA (057) , ARTRRA (058) , ARTRRA (059) , ARTRRA (060) 

• -5.10, -4.40, -3.90, -3.40, -3.10, 

• -2.70, -2.60, -2.60, -2.70, -2.40/ 

DATA ARTRRA (061) .ARTRRA (062) .ARTRRA (063) .ARTRRA (064) .ARTRRA (065) 

• ARTRRA (066) .ARTRRA (067) .ARTRRA (068) , ARTRRA (069) .ARTRRA (070) 

• -2.90, -3.10, -3.60, -4.10. -4.90, 

• -5.40, -6.40, -7.40, -8.20, -9.90/ 

DATA ARTRRA (71)/-11. 10/ 

THIS 50BR00TXRB HBTORRJ* A 6AIH VALUE FOR AR AR6LE 

AROLE IS THE DIPFHRBRCE OF THE ARTBRRA ARGLE AMD THE PITCH 
ARTARC IS THE A ^ TERRA AROLE WITH RESPECT TO THE AIRCRAFT 
PITCH IS THE r. ? .*CBAPr PITCH ARGLB 
ALL ARGLES AR5 EXPECTED TO BE EXPRESSED IR PADIARS 

THE RARGE OF VALVES FOB ARGLE IS: -70.9 TO 70.9 DEGREES 

IF AR ARGLE FALLS OOT OF THIS RARGE, THE CLOSEST GAIN 
VALUE MILL BE RETORMBD 

ARGLE - ARTARG-PXTCH 
IRDEX* IFXX ( (AHGLB*57. 3*70. ) /2. ) ♦ 1 
IP(XBDEX.LE.O) GO TO 30 
IP(IRDBX. GE. 71) GO TO 10 
ARGLE* FLOAT (IRDBX-1) *2. -70. 

GAIH*ARTRRA (IRDEX) ♦ (ARTRRA (IRDEX* 1) -ARTRRA (IRDEX) ) *ABS 
S ( (ARGLE*57. 3- ARGLE) /2.0) 

RBTORR 

10 GAIR*ARTRRA(71) 
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COMPARISON OF SCATTEROMETER RESOLUTION 
CELL AREA APPROXIMATIONS 
by 

Billy V. Clark and Richard W. Newton 


1.0 GENERAL 

For some time there have been two approaches to evaluating area, 

A, used In the scattering coefficient equation 

a m 

= X* A Kt G-rGt 

The approach used by NASA/ESD was to simply calculate the rectangular 
area along the ground track line. Figure 1. To account for roll and drift 
angles this area was divided by cosine roll times cosine drift. 

The approach used at TAMU is to compute the actual cell length between 
doppler contours, BB' in Figure 2, accounting for actual roll and drift 
values. The width of the cell is approximated very closely by the expres- 
sion for W* in Figure 2. The product, BB'xW, is the approximation to 
the area that Is used by TAMU In scatterometer data reduction. 

Since there have been some questions about the relative accuracy of 
the two area calculations, more detailed analysis of the area boundaries 
and more refined calculations of the area have been performed at TAMU 
Several types of calculations were made: 


BWA = ACTUAL BANDWIDTH 


(BWA)XH cos^ 
2 V cos^ 0 L 


FDL FDU 



2H tan 

W' * cosGcos^cosy 
Area - LW* 

Figure 2. Equations used by NASA/ESD to compute area 
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Figure 2. Equations used by TAMU to compute area. 


a. Area of trapezium (no parallel sides) 

b. Area of N parallelograms 

c. Integration by Simpson’s Rule 

d. Quadratic Formula to calculate BB*. 

The results of these different approaches to calculating the ground cell 
area Indicate In all cases that the technique currently used at TAMU Is 
a better over-all area estimator than the rectangular approximation. In 
the course of this analysis the mathmatlcal development of the expres- 
sion for y, page 9 of Reference [1], was reviewed. It was found that 
If the sign used before the square-root symbol Is positive, better agree- 
ment was achieved with the more accurate area values found by numerically 
Integrating over the cell. The form of the equation that should be used 
follows: 




4 
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2.0 AREA ESTIMATE BY TRAPEZIUM 


This approach makes only one assumption In computing the area; l.e., 
the lines connecting the four corner coordinates are straight* see Figure 3. 
Examination of the actual coordinate values at A, B, C, and D for viewing 
angles between -4° and -52° confirm that this Is a reasonable assumption. 

The calculation technique Is to use the equations for the four lines 
labelled XU, XL, FL and FU, Figure 2, to derive the actual coordinates at 
A, B, C, and D. The value for X at any point along XU can be shown to be 
given by: 

* (n-fi < 2 > 

where 

PI * COS(pCOS\p 

P2 • s1n$cos<|i 

M * (tan B/2)/s1ne' L 

8 = antenna bemwldth 

referenced to the aircraft 

0', » viewing angle through the 

L antenna coordinate system 

H * aircraft altitude 

4 • drift angle 

i|> ® roll angle 


H* 




Similarly, along the curve XL, the value of X Is given by 


„ _ > {rt-tfcMfh/tsmr (3) 


And, of course, along the doppler contours the value of X Is given by 



/»»*& J 


(3.*) 


where * viewing angle referenced to the ground track, Y-axis. By 
Iterating the value of Y, beginning at YL, the corner coordinates are 
found. 


Next, with the assumption 

that the connecting lines are straight, the 

area Is found by 





A 

» (A1 + 

M) 

- (A3 

♦ A4 + A5) 

(4) 

where 
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CM 

X 

(6) 

A3 

* (v c - 

V 

tX C 
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(7) 

M 

• <v 

V 

<V 

X a )/2 

(8) 

AS 
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< X 0' 

X A> 
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3.0 AREA ESTIMATION BY STRIP INTEGRATION 


This technique uses the sum of the area In N strips, each consisting 
of a selectively small parallelogram, to estimate the area. As Illustrated 
In Figure 4, the length of mini -rectangle Is the distance between the 
upper and lower frequency doppler lines. The length, DL, Is calculated at 
each step using the expressions 

>Li * ()0) 

(/- Kk. e *»*♦,) 

y u; - +[h<?'i'i(K u -o+K u .to*P-{l*') (n) 


'*** Kl* ( 2 % fS 

(12) 

<* = 

(13) 

and DL - (Yu. - Yl)/^P 

(14) 


The f4rst step was performed with ^ ° ty-6/2* and at each successive step 
with ^+1 8 ^ using * 0/500? 


111 



Figure 4. Area of N Parallelograms. 
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The value of DW^ was calculated by first evaluating range, 
the center of the line 01 and then using the expression 

J>W{ s frfacAfeea fi (l.ti 

The bc$i'> for ♦‘•e above expression Is shown In Figure 5. 



3>x/; xn*4to *•?• * A 

j>m/* %+*.4<r> m&Q. +4*tW&»(V; + a?)) 

Figure 5. Evaluation of EM 
UP 


The value of DU 1 was calculated by first evaluating range, 
the center of the line DL and then using the expression 

j>Wt & facAfe** fi (i. **»*M«*«3!) 

The basis for the above expression Is shown In Figure 5. 
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Figure 5. Evaluation of DW 
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4.0 AREA ESTIMATION FOR NUMERICAL INTEGRATION (SIMPSON'S RULE) 

This approach Is similar to the trapezium technique, but makes use 
of a more refined algebraic expression for areas beneath a curved line 
( 2 ). 

Using the equations (2), (3), and (3. a) above, and having already 
found the Intersection YL, YU, A, B, C, and 0 In Figure 3, repeated 
applications of Simpson's rule for numerical integration are made to 
find the areas: 

Al: beneath line FL from Y = Y(A) to Y = Y(B), 

A2: beneath line XU from Y * Y(B) to Y = Y(C), 

A3: beneath line FL from Y = YL to Y = Y(A) , 

A4: beneath line XL from Y = Y( A) to Y = Y(D) , 

A5: beneath line FU from Y * YU to Y * Y(D), 

A6: beneath line FU from Y = YU to Y = Y(C). 

For Al, A3, A5, and A6 a value of N - 8 was used. For A2 and A4 a value 
of N * 32 was used, where N specifies the number of times the Interval 
is halved between the limits of integration. The total cell area is then 
calculated as: 

A * Al + A2 + A5 - (A3 + A4 + A6) 


12 3 




5.0 QUADRATIC FORMULA APPLICATION 


Using the expression for Y In the quadratic form given at the top of 
page 9 t reference [1], the value of Y1 and Y2 were calculated. Next, 
the final expression on the same page was Implemented, using the value of 
U' specified on the same page. 


J0ar~ 


6.0 COMPARISON OF RESULTS 

The results of area calculation using the above techniques were 
compared to the results of the Rectangular and Doppler techniques. ‘These 
are sumnarlzed In Table 1. The results are graphically presented In 
Figures 6 and 7. The tables clearly show that the TAMU method has better 
general agreement with the Integrated values for area. 

Notice that In all cases where there Is a drift value the Rectangular 
area approximation departs rapidly from the other estimation techniques. 
Although the so called Doppler Contour approximation Is not considered 
perfect. It appears to be the best of the two "short cut" methods, Figures 
1 and 2, for handling excursion In drift. 
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TABLE 1 

Area Comparisons for Sever Calculation Tecnlques 

AT * Trapezium (4 - corners connected by straight lines) 
AP * Parallelogram (integrating of stripes) 

AQ * Quadratic Formula 

AR » Rectangular (ESD algorithm) 

AD s Doppler contour (Reference [1], page 9) 

Rol 1^ • 0 . 


Viewing 


Drift 
0, Deg 

Angle 
e L , Deg 

AT 

AP 

AQ 

AR 

AD 

AR— AD 41 nn 
“Id" 100 

0 

-5 

2104.8 


2100.2 

2100.2 

2100.2 

0.0 


-10 

2207.4 

2211.0 

2202.4 

2202.4 

2202.4 

0.0 


-20 

2222.1 

2211.35 

2209.6 

2209.6 

2209.6 

0.0 


-40 

3053.4 

3027.06 

3024.7 

3024.7 

3024.7 

0.0 

-4 

-5 

1982.0 


2105.5 

2105.3 

2105.5 

0.01 


-10 

2161.4 

2216.7 

2208.5 

2208.0 

2208.5 

-0.02 


-20 

2209.9 

2225.0 

2217.9 

2215.7 

2217.9 

-0.10 


-40 

3001.7 

2984.1 

2974.6 

2959.6 

2974.6 

-0.50 

8 

-5 

1907.1 


2121.5 

2121.0 

2121.5 

-0.02 


-10 

2141.6 

2251.7 

2217.3 

2215.4 

2217.3 

-0.09 


-20 

2232.0 

2267.0 

2243.0 

2234.2 

2243.1 

-0.40 

t 

-40 

3092.6 

3091.7 

3058.8 

2996.5 

3058.8 

-2.04 


If 


7.0 EVALUATION OF ROLL EXCURSION ON AREA VALUES 


To Illustrate the Influence of roll on the ground cell area a 
typical viewing angle, -45®, and a typical drift angle, 5®, were selected. 
Successive ground cell areas were calculated for successively larger roll 
angles, beginning at 0.0®, and remaining on the same doppler center fre- 
quency . Table 2 summarizes the results. Again, as was the case In drift 
excursions, the rectangular area technique does not adequately approximate 
changes In area for the off -ground-track situation. 

Finally, as can be seen from Table 2, the difference between the 
TAMU Doppler area and the Integrated area, by Simpson's rule. Is never 
greater than 0.5% for excursions In roll up to 6.7 degrees. 
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TABLE 2 

Area Calculation Comparisons 

« 

AD * Ooppler Area 
AR * Rectangular 

AQ • Quadratic Formula Application 
AS * Simpson's Rules Integration 


Roll, 4) 
Deg 

DRI FT, <J> 
Deg 

AS 

AQ 

AD 

AR 

*ioo 

0 . 

5. 

3132.9 

3132.9 

3133.9 

3074.4 

-1.9 

-1.2 

5. 

3157.7 

3141.9 

3153.9 

3075.1 

-2.5 

-2.5 

5. 

3184.1 

3166.3 

3178.4 

3077.2 

-3.2 

-4.9 

5. 

3241.7 

3228.8 

3241.2 

3085.8 

-4.8 

-6.7 

5. 

3289.1 

3287.7 

3300.4 

3095.8 

-6.2 


1* 
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R. W. Newton and B. V. Clark 


INTRODUCTION 

The Remote Sensing Center has been funded by The Earth 
Observations Division of NASA/JSC through Contract NAS9-1487S 
to evaluate the radar scatterometer processing techniques 
developed by TAMU to the one recently implemented by the 
NASA/JSC Engineering Systems Division. After several months 
of effort and many personal and verbal communications, the 
two processing techniques were demonstrated to be comparable. 
The only subsystems that were not satisfactorily compared 
were the area computational procedure. TAMU uses the algo- 
rithm documented by Schell in Technical Memorandum RSC-119 
[1], [2]. NASA/JSC used another algorithm. The largest 
discrepancy between these algorithms occurred at large roll 
and drift angles [3]. In order to insure that the output of 
each processing system (TAMU, NASA/JSC, hardware processors) 
were as consistant as possible, it was mutually agreed that 
the same area algorithm should be used by both TAMU and 
NASA/JSC. The algorithm that was agreed to by both TAMU and 
NASA/JSC is the one utilized by Schell (!]. It is the purpose 


of this memorandum to clarify the development and utilization 
of these equations. 


i 

l 


| 


i 


CELL LENGTH COMPUTATION 

The area algorithm approximates the true cell area by 
computing the length of the cell along the antenna center and 
multiplying this length by the width of the cell at the cell 
"center*’ defined by the filter center frequency. The cell 
length at each incident angle is computed by determining the 
points of intersection of the antenna center line with the 
isodoppler hyperbolas defined by the center frequency and 
bandwidth of the filter corresponding to each incident angle 
(Figure 1). The equation of the line going through the center 
of the antenna beam projected onto the surface can be deter- 
mined by its slope and intercept (Figure 2). The equation 
of this line is: 

X * (slopcU V(xci*.> 

From Figure 2 it can be seen that the slope is the tangent 
of the drift angle. For the geometry in Figure 2 

Si apC ~ 4 

Note that the drift, 4> , is considered to be positive when 
the nose of the aircraft is to the left of the aircraft 


- A - 

1 33 



cell length: B to B' 
cell width: C to C 

Approximate cell area: cell length x cell wid' 
Figure 1. Definition of cell length and cell width. 


s 


center of antenna beam 



Aircraft 

ground 

track 



<}> - drift angle (positive when aircraft nose is to left of 
ground track) 

x = (slope) y + x Q 

Figure 2. Geometry describing projection of the center of the 
antenna beam onto the surface. 


ground track. The x axis intercept, x Q , is given as 


H f 

*• " co* 4> * 

where the angle, 4* , is measured off of the z axis and orthogo- 
nal to the center line of the antenna. It is considered to be 
negative when the left wing goes low. The equation of the 
line defining the center of the antenna beam is now 


w wv - 

X = 4 . >3 - ^ 


( 1 ) 


The cell length is determined by computing the inter- 
section of this line with the isodoppler lines that define 
the cell end boundaries. The equation defining the isodopp- 
lers is determined from the doppler equation 



r 





where v - the velocity vector (in the direction of 
the y axis 

r - the range vector from the aircraft to the 

cell center (thus — — is a unit vector 

If I 

denoting the direction from the antenna 
to the cell center) . 


doppler frequency 


but from Figure 3 


V • ? * W l If \ u^a © r „ 

»\vllcl oa^Ce' t ^ 

V • c - \^| if 1 -bvw © / 

where 0 rv is the angle between the velocity vector and the 
range vector. Now 

ZlvllCl 

^ a N \f t 

Note that 

\F \ Sa/v\ 8' - - «j 

However, for small roll and drift angles 

^vv\ © X © 
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137 



making this assumption , then 


\F\ W*V 6 - 


and 

A \f \ 

Substituting for |r| 
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The factor k is constant for each incident angle. 
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The cell length can now be computed from (Figure 4) . 

cell length ■ - — (3) 

CO*. 4 

The y values are computed by solving for the intersection 
of the line defined by equation (1) and the hyperbolas 
defining the isodopplers, equation (2). Eliminating the 
variable x from these equations yields 
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a^ l + s.^ vc =t> 


'i 


- \o + i - 4ac? 




• $ ' 





CELL WIDTH 


The cell width can be conputed with the aid of Figure 5. 


VL 


W '* r tew 






The cell width projected onto the surface is 


2.V* tovA. 

VM ' 4 

Usxb COxf 


(4) 


AREA COMPUTATION 

The approximate cell area is computed by simply multi* 
plying equation (3) by equation (4) . 


Area ■ 



Area • 


ask.H' V. / 


(5) 
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Looking down the center line of the antenna. 



Figure 5. Cell width geometry. 





REDEFINITION OP ROLL ANGLE 


If it is determined that the roll angle is measured 
off the true nadir in the x z plane instead of orthogonal 
to the centerline of the antenna beam as assumed above, 
equation (5) can still be used to compute the area. However, 
the roll angle 9 used in the equations must be computed in 
the following manner: 

f 9 Tc*m * ( 


where ♦ is the angle in equations (1) through (S) above, 
and ♦ ' is the roll angle obtained from the NERDAS. 


COMPARISON TO NASA EQUATION 

The quadratic equation used to solve for »»nd is 


-tin 






where 



( 6 ) 


It can be shown that this equation is identical to the 
version that is contained in the NASA ESD program. The 
quadratic in the ESD equation is 




') i tts^©Cc K 1 4 -ll «-t*H <*4 

+ *&h *Q =o 

The constant k is obtained from 


$<S N 



Substituting this into equation (6) gives 

^ Wlfl - v* ^ ^ J a 

'i 1 'C'tSr)] - l* rtt5 -+ *• "] =o 

*"4^4 = ° 

yX -Q - '*’ H t *- lfe -l-coi , -4 f 0 "° 



This equation can be seen to be identical to the quadratic 
that BSD uses. 

CONCLUSION 

The equations developed in this document are identical 

i 

to those developed by Schell in Technical Memorandum RSC*119 
[11. A question has recently been raised by NASA BSD concern* 
ing the reference for the roll angle that is measured by the 
NERDAS. It has been shown that the equations developed in 
this document can be used for either roll angle reference 
with a very simple roll angle adjustment. 

The equations developed by Schell [1] are used in the 
real time 13.3 GHz radar processor and have always been used 
in the TAMU software processing programs. The only modifies* 
tion to the TAMU area computation that has ever been made is 
a change in the sign on the radical associated with the solu* 
tion to the quadratic equation. 
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